The quark-gluon vertex in Landau gauge QCD: 
Its role in dynamical chiral symmetry breaking and quark confinement 



Rcinhard Alkofer 1 , Christian S. Fischer 2 , Felipe J. Llanes-Estrada 3 and Kai Schwenzer 1 
Institut fur Physik der Karl-Franzens Universitdt Graz, 
A-8010 Graz, Austria. 
2 Institut fur Kernphysik der Technische Universitdt Darmstadt, 
D-64289 Darmstadt, Germany. 
3 Departamento de Fisica Teorica I de la Universidad Complutense, 
28040 Madrid, Spam 

(Dated: April 18, 2008) 

The infrared behavior of the quark-gluon vertex of quenched Landau gauge QCD is studied by 
analyzing its Dyson-Schwinger equation. Building on previously obtained results for Green functions 
in the Yang-Mills sector we analytically derive the existence of power-law infrared singularities for 
this vertex. We establish that dynamical chiral symmetry breaking leads to the self-consistent 
generation of components of the quark-gluon vertex forbidden when chiral symmetry is forced to 
stay in the Wigner-Weyl mode. In the latter case the running strong coupling assumes an infrared 
fixed point. If chiral symmetry is broken, either dynamically or explicitely, the running coupling 
is infrared divergent. Based on a truncation for the quark-gluon vertex Dyson-Schwinger equation 
which respects the analytically determined infrared behavior numerical results for the coupled system 
of the quark propagator and vertex Dyson-Schwinger equation are presented. The resulting quark 
mass function as well as the vertex function show only a very weak dependence on the current quark 
mass in the deep infrared. From this we infer by an analysis of the quark-quark scattering kernel a 
linearly rising quark potential with an almost mass independent string tension in the case of broken 
chiral symmetry. Enforcing chiral symmetry does lead to a Coulomb type potential. Therefore we 
conclude that chiral symmetry breaking and confinement are closely related. Furthermore we discuss 
aspects of confinement as the absence of long-range van-der-Waals forces and Casimir scaling. An 
examination of experimental data for quarkonia provides further evidence for the viability of the 
presented mechanism for quark confinement in the Landau gauge. 

PACS numbers: 11.10.St,11.15.-q,11.30.Rd,12.38.Aw 



I. INTRODUCTION 



Despite the progress seen in the last years the mechanism for confinement of quarks and gluons into colorless 
hadrons still lacks a detailed understanding. Based on the old idea that infrared divergences prevent the emission of 
colored states from color-singlet states, see e.g. refs. 0, different functional approaches have been employed to 
study the infrared behavior of QCD Green functions. It has turned out that Landau gauge is hereby the most suited 
one. Restricting oneself to the pure gauge sector it is even possible to make definite statements about the infrared 
behavior of all n-point functions @, Of. 

An important result of these investigations is the fact that the running coupling as determined from the Yang- 
Mills vertex functions shows an infrared fixed point 0, 0, 0, H| • This is obtained by balancing infrared suppressed and 
infrared divergent proper Green functions: Depending on the number of ghost and gluon legs, many of the one-particle 
irreducible Green's functions are indeed infrared singular, others are infrared vanishing. 

Hereby especially the infrared suppression of the Landau gauge gluon propagator immediately raises the question 
about the cause of the linearly rising quark potential when expressed in terms of Green functions. Additionally, 
dynamical breaking of chiral symmetry requires an effective strength of the quark interaction which, with an infrared 
suppressed gluon propagator, necessitates an enhanced quark-gluon vertex [9(. 

Therefore the quark-gluon vertex as the central link between the Yang-Mills and the quark sector of QCD is at 
the focus of many contemporary studies. It is also an important ingredient into the quark-antiquark interaction that 
is responsible for the formation and properties of bound states. On a perturbative level, the quark-gluon vertex has 
been studied in detail in arbitrary gauge and dimensions in [icj |. In this work, however, we are interested in the 
nonperturbative properties of this vertex. Lattice studies have granted some first insights into the relative importance 
of several components of the vertex at intermediate momenta [ill . [T3 |. In an attempt to understand these lattice 
studies some of us [l3|, [l4| and others [l5|, [l6|, [I?) have presented model studies based on resummations of one-loop 
results. As we will explain below such a treatment is not sufficient, consequently, a fully self-consistent treatment is 
required. 

An, at least at first sight, surprising result of such a self-consistent treatment is that dynamical chiral symmetry 
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breaking takes place in quite a non-trivial way also in the quark-gluon vertex: This vertex non-pcrturbatively assumes 
components of a Dirac scalar type in a similar self-consistent fashion as the quark propagator does. Especially, the 
assignment of a definite Lorentz structure is then not plausible. In view of the fact that the Lorentz nature of confine- 
ment, apparently indicated by phenomenological spin splittings, has worried two generations of QCD practitioners 
(and is still of interest 0, EOj HHj H3 ) this result appears, a posteriori, less surprising. 

One of our important findings to be reported in this paper is that these parts of the quark-gluon vertex have a 
considerable strength in the infrared. In particular infrared singularities arise not only in the uniform limit when 
all external momenta tend to zero, but also when only the gluon momentum vanishes. Interestingly, this soft-gluon 
singularity is, in contrast to those in the Yang-Mills sector j23[, as strong as the uniform divergence. This implies, of 
course, a complicated structure of the effective low-energy quark-quark interaction with the potential to permanently 
confine quarks. The remarks made above provide more than ample motivation for a concise study of the Landau 
gauge quark-gluon vertex. As will become evident below such an investigation reveals many unexpected qualitative 
features. 

This paper is organized as follows: In section II we introduce our notations and discuss the structure of the functional 
equations for the quark-gluon vertex. In section III we then develop an analytical infrared analysis for the vertex-DSE 
built upon our previous experience with the Yang-Mills sector of QCD. We discuss selfconsistcnt power law solutions 
of these equations in the cases of broken and restored chiral symmetry, investigate the corresponding running coupling 
and determine possible dependencies of the infrared solutions on the quark masses. In section IV we then present 
our numerical results thus confirming our analytical study. The infrared divergent behavior of the quark-gluon vertex 
found in our analysis implies a linear rising quark-antiquark potential extracted from the quark four-point function. 
This is demonstrated in section V. In the following section VI we deal with theoretical aspects of this confining 
potential like Casimir scaling and N-ality of the string tension. Phenomenological aspects are discussed in section 
VII, before we end with with a summary and conclusions. Many of the details of the corresponding calculations are 
relegated to various appendices. 

II. THE QUARK-GLUON VERTEX DYSON SCHWINGER EQUATION 

The fundamental dynamical clement in our improved analysis of the quark dynamics in quenched QCD is the 
quark-gluon vertex. In this section we discuss its basic definition, the functional equations it obeys as well as suitable 
approximation schemes. 

A. Basics 

The dressed quark-gluon vertex 
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is given by a combination of the twelve independent tensors that can be built from the two independent momenta p\ 
and P2 and the Dirac gamma matrices. Since the pioneering paper by Ball and Chiu [2~i| it is somewhat standard to 
divide this basis into a part transverse to the gluon momentum, Ti—x,.&, and a so called "longitudinal" part 1-h—x.a- 
The latter is not really longitudinal as it contains both longitudinal and transverse parts. The resulting basis is 
complete but not orthonormal and is explicitly given in eq. (|A lj) in the appendix. It contains the tree level tensor 
7 M that is a vector in Dirac space as well as scalar parts like i(pi + P2) 11 that are not chirally symmetric analog to a 
scalar mass term in the quark propagator. Introducing a set of scalar functions Ai, t; of the three scalar variables p\, 
v\i (j>2 ~ Pi) 2 , the vertex reads 

i=l i=l 

In the course of this paper we will give numerical results for the dressing functions Xi and r^. However, to demonstrate 
certain properties of the vertex at some places we resort to ansatze for these functions. To avoid confusion in the 
diagrammatics we represent the quark-gluon vertex by an open circle 
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whenever we use an ansatz instead of the fully calculated vertex given in eq. (flj). 

In addition we use the following diagrammatic expressions for the full propagators of QCD. The dressed quark 
propagator is given by 



* • 3 = ( i/ p* v (p 2 ) + a s (p 2 ))5 c ° lm (3) 

V 

with the vector part a v (jp 2 ) and the scalar part a s (p 2 ). These functions can be expressed in terms of the quark mass 
function M(p 2 ) and the quark wave function Zf(p 2 ) via 

a W) , n ,2, Z f (p 2 )M(p 2 ) 

ff,(p) -^iMV)' asip)= W+mW) (4) 

that are also expressed in this work in terms of the functions A{p 2 ) = 1/Z(p 2 ) and B(p 2 ) = Z(p 2 )M(p 2 ). The dressed 
gluon propagator is given by 

%TOBwromr =p»M^C l 5T t (5) 

p p v p z 

with the gluon dressing function Z(p 2 ) and the transverse projector P^ u {p) = (j)^ — ^-pr 1 )- Finally, the ghost 
propagator is given by 

a . b =-^r. (6) 

H p v P z 

It should be stressed that throughout this manuscript all propagators appearing in loop diagrams are to be taken as 
dressed, although we do not write the full circles in these cases so as not to clutter the diagrams. We also frequently 
refer to higher Green's function with more than three external legs. Here we generally denote one-particle irreducible 
(1PI) and connected functions by filled circles and ellipses, respectively, whereas connected or two-particle irreducible 
functions (2PI) are denoted by filled rectangles. 



B. Functional equations 

The Dyson-Schwingcr equations (DSEs) are given as an infinite tower of coupled nonlinear integral equations, sec 
[25l. l26l. I27I. 28| for reviews. Therefore most explicit computations require a specification of an approximation scheme. 
For such a scheme it is desirable to respect all symmetries of the underlying theory 1 . 

The simplest DSE in the quark sector is the one for the quark propagator given in fig. [T] which involves one bare 
and one dressed quark-gluon vertex. This equation has been studied in numerous works in the past. Of particular 
relevance for our work is ref. Q , where it has been shown that the vertex needs to be enhanced in the nonperturbative 
momentum region to provide for dynamical chiral symmetry breaking. 



1 Since the Lagrangian density with minimum derivative number and all the symmetries of QCD is QCD itself, any truncation of the 
equations of motion either 1) breaks at least one symmetry, or 2) is not derivable from a Lagrangian density, or 3) includes higher 
derivative terms. 
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FIG. 1: DSE for the quark propagator including one bare and one proper quark-gluon vertex. 

The quark-gluon vertex satisfies a full Dyson-Schwinger equation given in fig. El (see @ and references therein). 
This equation couples the proper quark-gluon vertex to higher connected four and five-point functions (the quark- 




FIG. 2: The DSE for the quark-gluon vertex where the r.h.s is written in terms of connected Green's functions. All internal 
propagators are fully dressed, also in all following diagrams. Note symmetry factors (1/2, 1/6) and signs have been omitted. 

gluon and quark-ghost scattering kernel, the two-quark-three gluon scattering kernel, and the quark-quark scattering 
kernel). In terms of one-particle irreducible (1PI) Green's functions the DSE, fig. [51 can be rewritten as given in 
fig. [3] These equations can also be directly derived in graphical form as outlined in [2!|. In both equations it is the 




FIG. 3: First version of the DSE for the quark-gluon vertex written in terms of proper Green's functions. 

external gluon leg that is connected to a bare vertex in the loop diagrams. Corresponding equations with external 
quark legs attached to the bare internal vertex exist. One possibility is given in fig. 3] Certainly, in the full theory 




FIG. 4: Second version of the DSE for the quark-gluon vertex written entirely in terms of proper Green's functions. 



both equations should give the same quark-gluon vertex. However, in a truncated theory, one of these two equations 
may be preferred over the other and more straightforward to solve. We deal with this point in appendix [Bl 
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One way to ensure that an approximation scheme respects at least those symmetries that are represented linearly, is 
to derive it from an action formalism. This holds for schemes based on n-particle irreducible (nPI) actions [30| . Such 
actions depend not only on the averaged field value but also on higher averaged n-point correlators thereby generating 
DSE-like equations via stationarity conditions. They are generalizations of the Cornwall- Jackiw-Tomboulis (CJT) 
action [3l|, [3^, [H, [H| which represents the order n = 2. Approximations based on the CJT action, which are also 
refered to as "^-derivable" , always involve fully dressed propagators whereas approximations based on an nPI-action 
involve dressed correlation functions with up to n legs. 

There are basically two ways to truncate nPI actions based either on loop or large- N expansions. The loop expansion 
of the 3PI action of QCD has been studied in [3(| up to 3-loop order. The result has the form of the standard CJT 
action which reads formally in Euclidean space 



where STr is a super trace over the generalized two point function T> for the supermultiplet <j) that includes the 
fields in the QCD action <fi = (A, c, c, q, q). The logarithmic term generates fully dressed 1-loop contributions whereas 
the 2PI-part of the action is represented graphically in Figs. 1-4 of ref. [3(|. Self-consistent equations for the 
Green's functions of the theory arc then obtained by variation with respect to the full propagators and vertices. The 
interesting result is that the truncated self-consistency equations are nearly identical to the DSE system truncated at 
the 3-point level when neglecting graphs with higher order vertices, but with the decisive difference that all vertices 
that enter the truncated self-consistency equations are dressed. We will later on confirm that this difference is not 



crucial when the 3PI equation is compared to the first version, fig. [3j of the quark-gluon vertex DSE. In section llTTl 
we will devise a truncation scheme for this DSE which reproduces the infrared behavior of both, the full DSE fig. [2] 
and the 3PI equation. However, things are slightly more complicated for the second version, fig. |4j of the quark-gluon 
vertex DSE. We discuss this point in more detail in appendix [5] where we argue that nevertheless all three equations 
lead to the same result in the infrared, in the presence of the "soft divergences" studied in subsection IIII El below. 

Finally we would like to mention that another important approach to 1PI Green's functions is given by the functional 
renormalization group equations (RGEs; for a detailed overview on this technique see p9l. l35l. l36l| ) . These have been 
successfully used in the past as an alternative technique to study the infrared behavior of Yang-Mills theory and 
implications for confinement 0, H3, [38| . The functional RGE for the quark-gluon vertex comprises of one-loop 
diagrams with all vertices dressed. Since this property is also captured by the nPI equations, we will not explicitly 
discuss these RGEs within this work. 



The quark-gluon vertex has been investigated in great detail in perturbation theory. Here Davydychev, Osland and 
Saks gave results for arbitrary gauge and dimension to one loop order and confirmed the validity of the corresponding 
Slavnov- Taylor identity However, as the momentum scale decreases and QCD becomes more strongly coupled, 
perturbation theory becomes inappropriate. One then has to resort to resummations and finally solve the vertex DSE 
self-consistently. As a first step it is useful to assess the importance of the different diagrams in the DSEs. This can 
be done within a semi-pcrturbative construction based on loop counting in the spirit of a particular type of skeleton 
expansion [39] in which all propagators are dressed but the internal vertices are still bare. This approximation can 
also be obtained from a two-loop approximation of the 2PI effective action. It allows us to assess the importance 
of the diagrams in fig. [3] on the basis of both, large N c counting and dynamical momentum dependences. In this 
approximation the DSE in fig. Q] is given by the diagrams of fig® which are conveniently labeled the "non-Abelian 
diagram" and the "Abelian diagram" .Whereas the latter diagram has a counterpart in the corresponding DSE for the 
fermion-photon vertex, the former has not. 




(7) 




FIG. 5: Equation of motion of the quark-gluon vertex as derived from a 3PI effective action. 



C. Comparison of Abelian and Non-Abelian diagram 
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FIG. 6: Semiperturbative approximation to the quark-gluon DSE. The first loop diagram is labeled the "non-Abelian diagram" 
and the second one "Abelian diagram". All propagators are dressed. 

These two diagrams and the corresponding semi-perturbative analysis is discussed in detail in Appendix[C]and allows 
to compare the relative strength of the Abelian and non-Abelian diagram numerically. From the color structure of 
the two diagrams one immediately notices that the Abelian diagram is suppressed by a factor of as compared to 
the non-Abelian one. Several authors [l5|, [l|| H3, 5l[ have recently employed truncations of the quark-gluon vertex 
based on the iV c -subleading Abelian one-loop diagram. This is necessary in the context of the Munczek-Nemirovsky 
model [42[ where one substitutes the gluon propagator by a Dirac delta function. Naturally, this operation is not 
possible in the non-Abelian diagram due to the product of two-delta functions. However those authors are aware 
that the Abelian construction is subleading by a factor of N% and has the wrong sign with respect to the dominant, 
non-Abelian construction. Therefore as an ad-hoc fix to alleviate the problem, they have changed the sign of the 
Abelian diagram and adapted the overall strength to the order of magnitude expected from the non-Abelian graph. 
The hope is then that the momentum dependence of both one-loop graphs is similar and therefore a calculation of 
the non-Abelian part can be sidestepped. 

We have checked this assumption; see Appendix [C] for details. Here we only discuss the result for a typical 




FIG. 7: Dressing factor Ai of the 7 M vertex in the semi-perturbative construction with dressed propagators and bare internal 
vertices. We show the total dressing function Ai = Zif + Xi A + X A (solid line, renormalized to 1 at fi = 13 GeV in the 
soft-gluon kinematic section) as well as the separate contributions Xi A and A^ from the non-Abelian and Abelian triangle 
loops. Finally, we show the non-Abelian contribution rescaled by an arbitrary factor to compare the momentum dependence. 



momentum configuration given in fig. [7] It can be clearly seen that differences remain even after arbitrary rescaling, 
so this rescaling factor must be dependent on the kinematics and adds to the systematic uncertainties. In general we 
find from our numerics that the N% suppression of the Abelian diagram w.r.t. to the non-Abelian one is even amplified 
by the dynamics. 

Another important difference between the Abelian and non-Abelian one-loop diagram (inherited by the semi- 
perturbative construction) is the quark mass dependence due to the quark propagators. The Abelian diagram contains 
two quark propagators, the non-Abelian only one. Therefore, upon increasing the quark mass, the diagram's leading 
pieces fall as 1/M 2 and 1/M respectively. This behavior can not be mocked by a constant factor unless this factor 
depends itself on the mass; but since different Dirac amplitudes depend differently on the quark mass, not all of them 
can be mocked simultaneously. 

Finally, we comment on the potential of this semi-perturbative construction to induce dynamical chiral symmetry 
breaking in the quark-DSE. To this end we inserted our semi-perturbatively calculated vertex (Abelian plus non- 
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Abclian diagram plus tree-level contribution) into the DSE for the quark propagator, fig. [TJ As a result we obtained 
only chirally symmetric solutions with B(p 2 ) = 0. This leads us to the conclusion that semi-perturbative trunca- 
tions of the vertex-DSE arc in general not suited to reproduce the patterns of dynamical chiral symmetry breaking. 
Consequently we now focus on a more selfconsistent approach. 

III. INFRARED ANALYSIS OF THE QUARK SECTOR 

After the general structure of the quark-gluon vertex DSE and its diagrammatic contributions have been studied 
in the last section via semi-perturbative methods, we will study its interplay in the full set of DSEs of QCD. This 
is done within a selfconsistent, purely analytical, infrared analysis of the coupled system of quark and vertex DSEs, 
which is later confirmed by our numerical solutions. 

A. Yang-Mills sector 

To begin, we shortly summarize previous results for the infrared behavior of the Green's function in Landau gauge 
Yang-Mills theory. We discuss the analytical behavior of these functions when all external momenta go to zero. 
Here, we choose all momenta proportional to each other and require the largest one to be much smaller than the 
typical scales of the theory, i.e. p 2 <C Aq CD . Here A qcD is generated dynamically in the gluon sector of the theory via 
dimensional transmutation. Then a self-consistent and unique power-law solution of the whole (untruncated) tower 
of DSEs is given by @, || 

r"<"V)~ ( P 2 ) [n - m) \ (8) 

where Y n ' m stands for the infrared leading dressing function of a Green's functions with n external ghost and antighost 
legs and m external gluon legs. The exponent k is known to be positive 0, EH but more restrictive limits have been 
suggested recently J23( ■ 

We stress that this IR behavior of gauge correlation functions is not influenced by the dynamics of massive quarks 
as will be confirmed subsequently. We note, however, that due to the possibility of strong divergences in the vertices 
this does not follow directly from their IR decoupling [ij [44[ . In the heavy quark limit even the value of the scale A QC d 
is determined completely by the gauge dynamics and obeys Aq CD <C M 2 . In contrast, in the chiral limit, where Aq CD 
is the unique scale of the system, there could in principle be a different conformal scaling solution in the Yang-Mills 
sector triggered by quark- loop insertions. This is discussed in more detail in (49j . Here we concentrate on quenched 
QCD, where ([5]) is certainly unaffected by the quark sector. 

Examples of the general solution ([8]) are given by the inverse ghost and gluon dressing functions, r 1,0 (p 2 ) = G~ 1 (p 2 ) 
and r°' 2 (p 2 ) = Z~ 1 (p 2 ), respectively They are related to the ghost and gluon propagators via eqs. (|5|6p . The 
corresponding power laws in the infrared are 

G{p 2 )~(p 2 y\ Z(p 2 )~(p 2 ) 2 \ (9) 

Since k is positive one obtains an infrared diverging ghost propagator, a behavior which is sufficient to ensure a well 
defined, i.e. unbroken, global color charge [5(J- 2 

Note that in Landau gauge an explicit value for n can be derived from the observation that the dressed ghost- 
gluon vertex remains (almost) bare in the infrared One then obtains K = (93 - Vl20T)/98 w 0.595 @, H EI, 

which implies that the gluon propagator vanishes in the infrared. A direct consequence of this behavior are positivity 
violations in the gluon propagator and therefore the confinement of transverse gluons 0, [53| . Another consequence 
of the solution ((5J) is the qualitative universality of the running coupling of Yang-Mills theory, see 3 and subsection 

riHGi 



2 We emphasize that the eq. (|8j solves the untruncated system of DSEs and the corresponding equations from the functional renormaliza- 
tion group. There is, however, an ongoing discussion, fostered by recent lattice results I451 I46H and also continuum investigations [47, 48] 
whether there is an alternative solution to the DSEs with an infrared finite or only weakly diverging ghost propagator and a massive 
gluon. This discussion involves subtile points of effects of Gribov copies and other artefacts on the lattice as well as in the continuum 
gauge theory. Pending further clarifications we will ignore this discussion for the moment and proceed by exploring the consequences of 
eq. ((S) for quenched QCD. We are confident that the results of our investigation, presented below also provides direct support for the 
validity of the solution eq. {8) . 
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Another important example of the solution ([8]) is the three-gluon vertex, which later on will serve as input for our 
analytical and numerical calculations. Restricting ourselves here to dressings of the tree-level structure r° (j>i,p 2 ,P3) 
(for a more detailed treatment see [H3]), this vertex goes like 

T^ Pl ,p 2 , P3 ) ~ (p 2 )- 3K Tl(PuP2,P3), (10) 

in all kinematical sections where there is only one external scale p 2 i.e. for p 2 ~ p\ ~ pi, ~ p 2 . 

As has been noted recently [23|, [54[ , there are additional kinematic divergences in the Yang-Mills 3-point functions 
that arise when only a single gluon momentum becomes small which we will refer to as the " soft-gluon limit" . The IR 
solution is parameterized by IR-exponents 6 for the leading correlation functions in the gauge sector. A lower index 
denotes the Greens function and the upper indices S u , 5 gl and S gh denote the uniform, soft-gluon and soft-ghost 
limits. The unique scaling solution is given in Table U It involves mild kinematic singularities in the three-gluon 



Sgh 


Sgl 


S 39 


S 3g 




S ah 








— K 


2k 





—3k, 


-4k 





1-2/t 


1-2/t 


1-2/t 



TABLE I: The IR exponents for the primitively divergent Green's functions of the unique IR scaling fixed point of Landau 
gauge Yang-Mills theory. 

vertex #j when only one external gluon becomes soft. 



B. Including quarks: basics 

Based on the infrared solution @ summarized in the last subsection we now extend our analysis to also include 
the quark sector of QCD. In the following we only consider quenched QCD and only briefly discuss possible effects of 
unquenching in section fVI Bl 

As a starting point we analyze the Dyson-Schwinger equation of the quark propagator, given diagrammatically in 
fig. [TJ Stripping the color factor, the quark propagator eq. ([3]) can alternatively be parameterized as 

ij> + M{p 2 ) 2 

g(p) = ^ + MV) /(P) ' { ) 

with the quark mass function M{p 2 ) and the wave function Zf(j> 2 ). In this notation the quark-DSE can be written as 

M( P 2 ) a f d*q M(q 2 )Z f (q 2 ) Z(k 2 ) 

- = Z 2 Z m m + g C F Z 1F tr / 77^-7^ 2 ; ,, 2 r 2 \ r *H fc ' Q) P t™ ,2 ( 12 ) 



Zj{p 2 ) ~ ^ J (2tt)4 q 2 + M 2 (q 2 ) vv ,HJ ^ k 2 

1 = Z 2+g 2 C F Z 1F tr ( ^Lj^ J Z { ( £\M k,q)P^ (13) 



Z f {p 2 ) ' ' ° r " J {2n) 4r q 2 + M 2 {q 2 ) " v ' p 2 k 2 

with k = p — q, the Casimir C F = (N 2 — 1)/(2N C ) and the renormalization factors Z 2 of the quark fields and Z\ F of 
the quark gluon vertex. The symbol tr denotes the Dirac-trace. 

Now consider the loop integral on the right hand side of the DSE. When chiral symmetry is broken dynamically we 
have M (p 2 ) ^ and the denominator of the internal quark propagator cannot be smaller than the generated quark 
mass M(p 2 ). As an effect the integral is dominated by momenta larger than the quark mass and the resulting quark 
mass on the left hand side of the equation freezes out to a constant in the infrared, in selfconsistent agreement with 
M(p 2 ) 7^ 0. This infrared screening mechanism of the quark, observed in a broad range of works so far (see Q and 
refs. therein), can only be overruled by an extremely singular quark-gluon vertex. We have not found a selfconsistent 
solution for the vertex-DSE that supports such an extreme case and therefore discard this possibility. In general this 
screening mechanism will also make the quark wave function constant in the infrared, as again observed repeatedly 
in the literature. 

The only other possibility we see is given by the 'scaling solution' Zf(p 2 )/AI(p 2 ) ~ (p 2 )~ K ~ 1 ^ 2 . Since this solution 
is refuted by lattice calculations of the quark propagator [5 51 ] we do not discuss it here, but refer the interested reader 
to [49[ for more details. 

We thus arrive at an IR finite solution of the quark propagator with 

M = AI(p 2 -> 0) and Z f = Z f (p 2 -> 0) (14) 
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FIG. 8: Dyson-Schwinger equation for the quark-gluon vertex. In the first line we show the full equation @). The second 
line shows important terms of the lowest order of a skeleton expansion, that already reveal the correct power counting in the 
infrared. All internal lines in the diagrams represent fully dressed propagators. Symmetry factors and signs have been omitted. 
The infrared counting analysis can be followed by referring to the second diagram in the right-hand-side of the second line, 
where one can abstract the ghost loop as if forcing a dressed three-gluon vertex. This is the same counting derived from the 
last diagram of figure [S] 



where we define the constant and non- vanishing limits by M, Zf and always refer to these limits in the following when 
no momentum arguments are given. The quark propagator can then be written as 



S(p) = 



ij> + M(p 2 ) 
p 2 + M 2 (p 2 ) 



Zf(p 2 ) 



M 2 



M 



(15) 



for momenta p 2 <C Aq CD . 

In the following we show that a selfconsistcnt solution of the quark-DSE and the DSE for the quark-gluon vertex can 
be found based on eq. (|15p . To this end we now apply the power counting method of [H, Q| to the Dyson-Schwinger 
equation for the full quark-gluon vertex, given diagrammatically in the first line of fig. [H In the second line the 
higher n-point functions of this equation have been expanded to lowest order in a skeleton expansion in terms of full 
propagators and vertices. The dressed quark-gluon vertex T M can be decomposed in a basis of twelve tensor structures 
[24| given explicitly in cq. (|A2[) in the appendix. 

For the internal loops of the skeleton expansion in fig. [8] there are two possibilities: 

(i) they are dominated by loop momenta of the order of the small external scale p 2 ; this then requires infrared 
singularities in the internal quark-gluon vertices, 

(ii) the same screening mechanism applies as in the quark-DSE; for a small external scale p 2 the loop diagram is 



then dominated by loop momenta in the mid-momentum region q 2 
out in the infrared. 



Aq CD and consequently the vertex freezes 



We emphasize that both possibilities may lead to selfconsistent solutions of the quark- and vertex-DSE, although 
numerically we only found solutions with the property (i). Whether this observation is truncation dependent or not 
needs to be investigated further in future work. Noting that the second possibility requires no further infrared analysis 
anyway we proceed to show that indeed selfconsistent solutions can be found when (i) is realized. 

In this scenario the dressing functions of internal Yang-Mills propagators and vertices can be approximated by 
the power laws given in eq. l[5]). Wc therefore use cq. for the scaling of the ghost and gluon propagators. The 
ghost-gluon vertex scales as a constant in the infrared [51(. The internal quark propagator lines are replaced by (|15p . 

We contracted the right hand side of the DSE with appropriate tensor structures to project out the various compo- 
nents of the vertex, c.f. eq. (|A2[) . Thus we eventually have to deal with scalar integrals only, which depend on powers 
of internal and external momenta. The evaluation of such diagrams has been described in [56| . Most important for 
our purpose is the fact that all powers of internal momenta transform to powers of external momenta after integration. 
This is also clear from simple dimensional consideration, provided there is only one external scale present. In our case 
this is ensured by the condition p 2 <C Aq CD < M 2 , which leaves only one scale p 2 in the deep infrared. The following 
analysis rests on the assumption that there is no second small scale present. In principle this assumption may not 
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be true for heavy quarks, since then one has a potential additional small scale Aq CD /M. We have investigated this 
possibility in our numerical calculations of the quark-gluon vertex and did not see any influence of such a quantity. 
The corresponding results of our quark mass study are presented in subsection IIVQ where we discuss this issue 
further. 



C. How not to include quarks 



Before we come to our actual solution we would like to point out the complication that the Dirac structure induces. 
Suppose we only consider a single form factor Ai (p 2 ) with a single IR exponent S qg for the quark gluon vertex corre- 
sponding to the tree- level tensor structure 7^. Due to the Dirac-algebra the integral then picks up only contributions 
from the vector part of the internal quark propagator. The resulting vertex-DSE (IR- leading diagram only), that can 
be followed from figure [5] and its caption, is then given by 

w„ v) _ ( N c a 2) f A Z(kl)Z{kl) Z f (q 2 ) 

A lM A l( ^)(f -Sl-^ + i ^) (16) 

with the momenta k\ = p\ — q, = q — P2 and the dressing T 0,3 of the three gluon vertex. Here we used tree-level 
tensor structure in the three-gluon vertex without loss of generality. After integration, all powers of internal momenta 
are converted into powers of external momenta for dimensional reasons (potential conversions into powers of the quark 
mass are subleading in the infrared and are therefore discarded here; we further justify this treatment a posteriori in 
our numerical section TlV C|) . 

Power counting in terms of the external momentum p 2 then results in 



(p ) qg ~ p A X 



P 2 



x ^ x (p 2 )- 3K * [MP 2 )]* * P 2 (17) 



M 2 



~ ^/ ,, 2 (18) 

Thus a selfconsistent power-law solution of this equation would entail that 

S qg = -1-K, (19) 

i.e. a singular quark-gluon vertex in the infrared. This solution, however, only serves illustrational purposes, since 
it does not survive when the full Dirac structure of the vertex included in the other tensor structures is considered. 
This is shown in the next subsection. 



D. Including quarks: the correct way 



To see how the Dirac structure of the vertex can introduce additional complications we note that the tensor 
structures of the vertex group into two different sets of Dirac-scalar and vector parts, as detailed in appendix I A II We 
have to take into account that within the DSEs the two parts receive qualitatively different contributions from the 
various loop corrections and introduce distinct IR exponents S qgv and S qgs for the vector (qgv) and scalar (qgs) Dirac 
structures. One then arrives at a coupled set of equations with various contributions: 

( ( 2\l + K+2S qgv ( 2\l + K+2S qgs ( 2\l/2+K+S qgv +S q 



M 2 1 J M 2 '~ J M 

( I 2\1/2+k+2(5 m „ I 2\1/2+k+25„ s 12 \ 1+K+S qgv +S qgs ~| 

(P 2 ) S - - — [ Zf (P J M - , Z, (P) M , Z f 1 (20 ) 

where the 'max' indicates that the most singular term in the infrared is the leading one on the right hand sides of 
these equations. This equation can be solved by a range of solutions given by 



Sqgv S [-1/2- K,-«], 5 qgs = -1/2 -K 



(21) 
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As is demonstrated explicitly in appendix [Dl this infrared solution for 8 qgv and S qgs does not change when also the 
tensor components of the quark-gluon vertex are taken into account. The above encountered behavior eq. (|19p is 
therefore an artefact due to the ommission of the scalar Dirac components. 

It is furthermore instructive to look closer into the mechanism that generates the solution (|2ip and recall the 
properties of the different tensor structures with respect to chiral symmetry. The corresponding transformations act 
on fermion fields via e lai tp with a labeling the transformation. A chirally symmetric quark propagator or quark-gluon 
vertex would therefore correspond to 

{ 7 5 ,r^} = and { 7 5 ,S*} = 0. (22) 

These relations are satisfied by all tensor structures with an odd number of gamma matrices. Those parts with an even 
number of gamma matrices can only be present when chiral symmetry is violated. From eqs. (|20[) we find that the 
driving pieces for the selfconsistcnt solution (|2"Tj) of the vertex DSE are the ones with broken chiral symmetry. If these 
structures are left out in both the quark propagator and the vertex, one finds a different (infrared weaker) selfconsistcnt 
solution, as discussed in subsection IIII Fl below 3 . The fact that the presence or absence of chiral symmetry breaking 
terms generates different solutions of the vertex-DSE seems a remarkable and novel result to us. It underlines that 
dynamical chiral symmetry breaking is not restricted to the quark propagator alone but happens (consistently) in 
every Green's function with fcrmionic content. Since the underlying mechanism is not well appreciated in the literature 
we found it useful to illustrate it with a simple mechanical model, discussed in appendix [El 

Before we discuss some aspects of the solution eq. ([2"Tjl , it is important to check that it persists to all orders in the 
skeleton expansion used in fig. [5] Higher order terms in this expansion can be generated by inserting diagrammatical 
pieces like 




(...) 



into lower order diagrams (all propagators are fully dressed). In ref. Q some of us showed that insertions involving 
pure Yang-Mills propagators and vertices do not change the overall infrared divergence of a given diagram. This 
is also true for the three pieces involving quark lines: the first one introduces two quark-gluon vertices, two quark 
propagators and one gluon propagator plus a new integration into a given diagram. This amounts to an additional 
power of 

/„2^-1-2k+2k-1+2 7c 

z ' J ^f' — =§<? 2 >° < 23) 

in the external momentum (p 2 ). A similar result is obtained for the other pieces. Thus diagrammatic insertions of 
quark lines do not change the infrared behavior of a given diagram but introduce additional powers of Zf /M 2 . This 
implies that the solution eqs. (f2"Tj) is valid to all orders in the skeleton expansion. It is therefore also an infrared 
solution of the full vertex-DSE. 

Note that similar to the DSEs in the Yang-Mills sector it is the diagram containing the ghost triangle loop that 
dominates the right hand side of the equation in fig. [H] Thus the infrared divergent behavior of the quark-gluon vertex 
is an induced effect from singularities in the pure Yang-Mills theory. 

Finally, we emphasize that the ghost triangle loop in the infrared leading diagram is also the leading infrared 
contribution to the fully dressed thrce-gluon vertex [3|. This offers an interesting opportunity for an approximation 
which enables our numerical analysis in section HVi Since we are mainly interested in the infrared behavior of the 
quark-gluon vertex we may replace the ghost triangle sub-diagram by the dressed three-gluon vertex as shown in 
fig. [9] The resulting diagram greatly resembles the "non-Abelian" diagram already discussed above in our semi- 
pcrturbative analysis, but with all bare vertices replaced by dressed ones. It is interesting to note that such a 
completely dressed non-Abelian diagram is also present in the equations of motion for the quark-gluon vertex from 
the functional rcnormalization group (FRG) and nPI formulations (with n > 2). This suggests that our infrared 
solution eqs. (|2"Tj) also solves the corresponding FRG and nPI equations self-consistently. We have checked that this 
is indeed the case. 



3 Note in this respect that the (inconsistent) behavior eq. I I19H has been obtained from the chiral symmetry breaking part of the propagator 
while the vertex had only the symmetry preserving term 7^. 
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FIG. 9: The infrared leading diagram of the quark-gluon vertex DSE is shown in the first line. In the second line we show an 
approximation of the full quark-gluon vertex DSE that preserves all infrared features of the equation, while also capturing the 
leading l/iV c -parts of the ultraviolet properties. 

E. Soft-gluon singularity in the quark-gluon vertex 

In addition to the main infrared singularity that appears when all scales in a given Green function are sent to zero, 
there can be kinematical singularities that appear in specific kinematic sections. The present counting rules relied 
on the implicit assumption that all external momenta pi scale as some common scaling variable p. However, in case 
any of the momenta vanishes identically this is not fulfilled. In the parameterization with a common scaling variable 
p and dependencies on the other momenta given in terms of dimcnsionlcss momentum ratios Pi/p € (0, oo) or angle 
cosines between momentum vectors Xi £ (—1,1), the IR expression reads 

F{ Pl /p lXl ){p 2 ) a , i = l...n-l 

and these kinematic singularities appear as singularities of the function F . We have identified such a kinematic 
singularity for the quark-gluon vertex in the soft-gluon kinematic section where the gluon momentum p% = p2 — pi = 0, 
whereas p\ and p\ may be large. In this configuration the singularity is even self-consistently enhanced and thereby is 
likely to constitute a stable solution of the DSE system. This singularity directly extends the uniform IR singularity 
discussed above and reported previously in (5?| to this larger kinematic section. 

The mechanism behind this singularity is based on the strong (p 2 ) K singularity in the three-gluon vertex when 
all its momenta vanish simultaneously. To see this consider the infrared leading diagram @ of the quark-gluon vertex 
DSE. 




We choose a momentum configuration, where the external momenta p\ and p\ may be large, but the external gluon 
momentum p^ = e is small. This allows to route the loop momentum <5 to expose the relevant phase-space where both 
internal gluon momenta, 8 and 5 + e, arc small. The incoming quark momentum is fixed at p\ and the outgoing is 
P2 = Pi + e: Now, if there are soft-gluon kinematic singularities in the quark-gluon vertex, these also appear in the 
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FIG. 10: Support of uniform versus soft IR divergences in the quark-gluon vertex. Left: if the quark-gluon vertex is expressed 
in terms of the three Lorentz invariants p\, p%, p 2 g = (pi — P2) 2 , the uniform divergence corresponds to the point where all 
momenta simultaneously vanish, whereas the soft divergence lies on a plane in this three-dimensional space, since only the 
gluon momentum is required to vanish. This shows the larger phase-space (and therefore larger contribution to loop integrals) 
associated with soft divergences. Right: if one uses instead of p g a Euclidean spherical angle between pi and P2, whose cosine 
x £ (—1,1) is on the OZ-axis, the uniform divergence corresponds precisely to this segment. However the soft divergence is 
now also a segment, at fixed x = 1, running for pi = pi- Thus a grid optimized to capture only the uniform-divergence in a 
numerical code will have difficulties to capture the soft-divergence too. 



internal quark-gluon vertices at this particular kinematical configuration and if these singularities are strong enough 
loop momenta around 6 will dominate the integral. 

Then, leaving out finite factors, self-consistency of the soft-gluon divergence demands that the corresponding expo- 
nent 5f g fulfills 

(e 2 ) 5 - oc J dH OS 2 ) 2 *" 1 {{6 + e) 2 ) 2 *- 1 (S 2 ) s ™ ((5 + e) 2 ) 5 - (S 2 + e 2 + (5 + e) 2 )- 3 * T° 3g (24) 

where (S 2 + e 2 + (5 + e) 2 )~ 3K stems from the dressing of the three gluon vertex. Note that the bare three-gluon vertex 
r® g also contributes a power of momentum (e 2 ) 1 / 2 after integration. 

Now, if e — > 0, the integral is inconsistently divergent in 5^0. Therefore, e needs to be kept finite, in which case 
the integral is regulated to the RHS of 

(e 2 ) 5 t oc ( e 2 )-+ 2 ^+V2 (25 ) 

and therefore it must be that 

5f g = -1/2 -k (26) 

which indeed corresponds to a soft singularity with an equally large exponent in the vertex. 

Below, we check explicitly that the counting of soft infrared divergencies is not affected by the remaining diagrams 
in the quark-gluon vertex DSE. In fact the situation is similar to the kinematic section where all external momenta 
go to zero: In a Dyson-Schwingcr formulation the non-Abelian diagram is the infrared leading one with all other 
diagrams being subleading. In an nPI formulation or within the FRG framework, where all vertices are dressed in 
general, we encounter a uniform scaling of the whole equation, i.e. all diagrams have the same infrared divergencies. 

The importance of the self-consistently enhanced soft-gluon divergence comes from its role for the confinement 
mechanism of quarks discussed below. A mechanism based purely on the uniform IR divergence as discussed above 
and proposed in [57j would require the external quarks to have vanishing Euclidean momenta - corresponding to 
light-like momenta in Minkowski space to trigger the divergence. This is only possible for chiral quarks, whereas any 
mass would act as a cutoff that shields the divergence. As known from the structure functions measured in deep 
inelastic scattering experiments [58], the quarks in hadrons have a broad momentum distribution which even vanishes 
at zero momentum for valence quarks. Similarly, constituent quarks inside charmonia or bottomonia have momenta 
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of the order of or larger than Aq C d (typically of order a s MQ). The soft-gluon divergence is realized independent of the 
quark kinematics and requires only the gluon momentum to become small. Therefore this type of singularities has 
the potential to provide a static linear potential between heavy quarks since the momentum of an exchanged gluon is 
always small as far as the quarks are sufficiently spatially separated. This will be discussed in more detail in sections 
Hand EH 

Let us now analyze the IR power counting for the full system of DSEs for the quark sector of quenched QCD. In 
contrast to the propagator which has only one external scale and is parameterized by a single IR exponent S q , the 
quark-gluon vertex has two independent external momenta. There are three different limits that can be IR sensitive, 
namely the uniform limit when all momenta vanish and the soft-gluon limit respectively the soft-quark limit when 
only one corresponding momentum vanishes. The corresponding IR exponents are distinguished by upper indices 
5™g, S^g and , respectively. Including all these possibilities in our analysis we can show that the uniform and soft 
divergences studied before indeed provide a solution of the full non-perturbative system of DSEs (for simplicity we 
display results only for the one-loop skeleton expansion, but the generalization to an arbitrary order should now be 
evident, and a skeleton-free proof of the dimensional counting along the lines of possible). 

In a case where scales of different order are present the IR scaling law can in principle involve all of them p3l | . Since 
it results from the conformal pure gauge sector, the part in the scaling law that involves powers of the anomalous 
scaling exponent k always involves the characteristic scale A qC D- When there are soft s and hard h external momenta 4 , 
the scaling law can also involve ratios of them and finally the quark propagators introduce the quark mass M. We 
write all scaling laws in terms of the soft external momentum over Aq CD and skip any finite ratios of the different hard 
scales. From the point of view of power counting the integral can be dominated by modes in the regions around all the 
external scales. Therefore, it has to be separated into two parts, namely an 'infrared' part where the loop momenta 
are of the order of the soft external momentum s <C Aq C D that scales to zero and an 'ultraviolet' part where the loop 
momenta are either of the order of the quark mass M > Aq CD or of the order of the hard momentum h < Aq C d that 
stays finite in the limit. 

The explicit power counting analysis is presented in appendix IFl and yields the following limits: 

HmS(p) ~ (p 2 ) 5 ' 
p— >o 

lim r„(p,«,A) ~ (k 2 +p 2 + q 2 f™ Tl(p,q,k) 

limr,(p,^) ~ (p 2 ) 5 ™ T°(p,q,k) 

limiyp^fc) ~ tfyh T°(p,q,k) 

g — >0 ^ 

\im T^p 7 q,k) ~ (fc 2 ) 5 « r°(p,g,fc) 

where k is the gluon momentum and p, q are the quark momenta. Here is the fully dressed quark-gluon vertex 
and T° are the undressed tensor structures given in eq. (|A2|) . Eq. ([27]) is the central result of our infrared analysis. It 
shows indeed that the strong dynamics in the quark sector can self-consistcntly induce both a uniform and a soft-gluon 
singularity in the quark-gluon vertex, that are of equal strength and which are fully compatible with each other. 

F. Restoring chiral symmetry 

The infrared power law solution found in the last subsection is valid in the chirally broken phase of quenched QCD. 
In the following we will study the chirally symmetric case, i.e. we will omit all those Dirac structures in the propagator 
and quark-gluon vertex that break chiral symmetry. We still restrict ourselves to the case of quenched QCD, i.e. we 
do not take into account the backreaction of the chiral quarks onto the Yang-Mills sector. The solution ((Sj) is then still 
valid. In general, however, unquenching effects are certainly important for chiral quarks and can lead to a different 
solution in the gauge sector. This possibility will be discussed in (49j . 



with Sq = 



With % = 1 1-2, 



with 8 q qg = 



with 



99 



l-2« 



(27) 



4 Note one could also speak perhaps more appropriately of "ultrasoft" and "soft" momenta but this might lead to confusion with NRQCD 
naming. 
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With chirally symmetric Dirac structures only, the quark propagator and quark-gluon vertex are given by 

SxiP) = tZf(p 2 ), T^ P )= t glj2Hp 2 )Lf+ £ ntf)T?\ (28) 

" \i=l,2 1=2,3,6,8 J 

in the presence of only one external scale p 2 ; the chirally symmetric vertex structures in the general case are given 
explicitly in the numerical section IIV El 

For the infrared analysis we project the vertex-DSE on the remaining tensor components. We then arrive at the 
reduced system for the IR exponents 5 qx and S qgx in the chiral limit 

( p 2)-^+i/2 „ max{(/) 1 / 2 ,(p 2 ) 2K+5 ^+^+V2j (29) 

(p 2 )<W „ ma^{(p 2 ) Q ,{p 2 ) K+2S ^ +s ^} (30) 

where Zf(p 2 ) ~ (p 2 ) Sqx denotes a potential power law for the quark wave function in the infrared. Assuming for the 
moment that Z(p 2 = 0) = Zf stays constant in the infrared, similar to the chirally broken case discussed above, we 
then find the solution 

S gx = °! S ggx = ~ K ' ( 31 ) 

i.e. an infrared finite quark wave functions and an infrared divergent vertex. Note, however, that the vertex is 
considerably less divergent than in the case of broken chiral symmetry above. 

The behavior eq. (|3~lj) is also realized in our numerical solutions, as shown in subsection IIV Bl We have also 
investigated the chiral quark-DSE for a range of possible behaviors of the quark-gluon vertex and found that this 
equation is always controlled by intermediate momenta. In the infrared this naturally leads to a constant Zf(p 2 = 
0) = Zt in agreement with our assumption from above. 

In addition to the solution presented above, a pure power counting analysis cannot exclude the possibility that both 
the propagator and the vertex scale trivially in the IR limit 



*™ = 0; <W = o. (32) 

This possible alternative fixed point is interesting since it would yield an IR vanishing running coupling according 
to the discussion in the next subsection. A similar behavior has been found previously for the scaling regime of the 
ungapped magnetic sector at high density [6(| ■ We defer a complete analysis of the fixed point structure for chiral 
quarks to a future work [49[ . 



G. The running coupling 

We now wish to investigate the running coupling a qg that can be defined from the quark-gluon vertex. To this end 
we first recall the corresponding definitions for the couplings from the Yang-Mills sector. Here, renormalization group 
invariant couplings have been defined from either of the primitively divergent vertices of Yang- Mills-theory, i.e. from 
the ghost-gluon vertex (gg) 7 the three-gluon vertex (3g) or the four-gluon vertex (4g) via 

%> 2 ) - ^-G 2 (p 2 )Z{p 2 ) p2 ~° const/iV c , (33) 

«3> 2 ) = |V<V)] 2 ^V) P -° const/iV c , (34) 

«4 9 (P 2 ) = f-r°<V)ZV) const/iV c . (35) 

Using the DSE-solution |8]) it is easy to see that all three couplings approach a fixed point in the infrared 0|. This 
fixed point can be explicitly calculated for the coupling eq. (|33|) . Employing a bare ghost-gluon vertex one obtains 
a gg (0) w 8.92/iV c 0]. A first step in the evaluation of the coupling of the four gluon vertex has been taken in [oil ], 
where the value a4 g (0) « 0.0083/A^ c has been found. 

A nonperturbative and renormalization group invariant definition of the corresponding coupling from the quark 
gluon vertex is given by 

a qg (p 2 ) = fx 2 1 (p 2 )Z 2 (p 2 )Z(p 2 ) p2 ~° 1, (36) 

47T J p 



16 



where Ai dresses the 7 M -part of the vertex. With the exponents S qgv = 5 qgs = — 1/2— k, as obtained in our numerical 
solutions, we find by power counting that the coupling is proportional to 1/p 2 . Thus, contrary to the couplings from 
the Yang-Mills vertices, we find this coupling to be singular in the infrared, which corresponds to the old notion of 
infrared slavery p], [6^, HH, that has been deemphasized for a decade or two. 

If, however, quenched QCD is forced to the chirally symmetric phase as done in the previous subsection, we recover 
the fixed point behavior also in the coupling from the quark-gluon vertex. Then the weaker infrared singularity (]31|) 
of the vertex leads to 

a qgx {p 2 ) = ^>?i{p 2 )Z 2 f { P 2 )Z{p>) p2 ~° const. (37) 

We thus uncovered a close connection between dynamical chiral symmetry breaking and two different possible behav- 
iors of the coupling from the quark-gluon vertex. One obtains infrared slavery when chiral symmetry is broken and a 
fixed point when chiral symmetry is restored. We will see later on in section [V] that this connection also extends to the 
quark-antiquark potential and thus suggests a novel mechanism linking confinement and dynamical chiral symmetry 
breaking in quenched QCD. 



H. Heavy-quark mass dependence of the quark-gluon vertex 

In the previous subsections we discussed in some detail the emergence of selfconsistent infrared power laws in the 
quark-gluon vertex. From our analytical infrared analysis we derived an infrared behavior of (p 2 )~ 1//2 ~ K in case of 
only one external scale p that goes to zero and a soft-gluon infrared divergence of (j>g)~ 1//2 ~ K in case only the external 
gluon momentum p g goes to zero. We now investigate the behavior of the corresponding coefficients and in particular 
their dependence on the quark mass. 



1. Naive counting: M(0) ~ m 

To this end we start with the vertex constructed in one-loop perturbation theory. Then the internal vertices are not 
dressed and only feature a color charge g that is universal and independent of the quark mass, with scale- running effects 
not yet manifest. In first-order perturbation theory, all dependence on the mass comes from the quark propagator, 
that can be written as 

a_ 1 ^+M _1_ Jt_ 

jk - M k 2 + M 2 m^oo M M 2 ' 1 ' 

Therefore, in the infinite mass limit, the contribution with an even number of Dirac 7 matrices (the scalar piece) is 
suppressed by one power of the heavy quark mass, the contribution with an odd number (the vector piece), by two 
powers. 

The behavior of the quark-gluon vertex in perturbation theory, given by fig. [S] where the dressed vertices are 
replaced by bare ones, is then easy to obtain. The vector piece is dominant thanks to the bare (mass independent) 
diagram, and subsequently takes 1/M 2 corrections from both the Dirac structure (7 fc-f)/M 2 of the non-Abelian 
(second to last) diagram and the Dirac structure (jMjMj)/M of the Abelian diagram. Likewise the scalar piece 
has no zero order term, and takes 1/M corrections from the structure (■jM'y)/M 2 of the non-Abelian diagram at 
one- loop. For brevity we will from now on use the symbol V\ for the vector part of the vertex (odd number of gamma 
matrices) and S\ for the scalar part (even number of gamma matrices), where the subscript indicates the number of 
gluon legs. Consequently we will use in the following Vb to denote the vector part of the quark propagator and So for 
its scalar part, i.e. we have 

Vo^l/M 2 , S ~ 1/M. (39) 

We will read Feynman diagrams following the spin line from left to right. In this notation we can state for the mass 
dependence at one-loop order 

Vi ~ 0(M°) + . 9 2 0(1/M 2 ) , Si ~ g 2 0(l/M) , (40) 

where the leading correction is in the scalar piece and arises from the ViSoVi-combination in the non-Abelian diagram. 

Notice that the Abelian diagram also gives sublcading corrections through the ViVoViVoVi-combination and the 
ViiSoViVbVi -combination. In terms of mass dependence these are further suppressed compared to the leading terms 
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(|40p and can be ignored. Thus the conclusion would be that in perturbation theory, heavy quarks increasingly see a 
Coulombic potential from perturbative one-gluon-exchange. 

We will now proceed beyond perturbation theory and discuss the leading N c , leading skeleton, self-consistent 
equation for the vertex as given in fig. [SJ Let us first consider the Dirac amplitudes with an odd number of gamma 
matrices. The bare term is of order 0(M°) whereas the one-loop skeleton graph on the right gives three types of 
contributions, namely VfVbVi, SiSqVi and SiVqSi. For the even amplitudes, as Si, the bare term does not contribute 
and the one-loop skeleton graph yields structures S±VoVi, V\SoVi and SxSqSx- 

The entire system is still consistent with the perturbative result, that is, with vertex counting 

Vi ~ M°, S*i - M~ l . (41) 

But now in addition self-consistent solutions arise with 

Vi ~ M°, Si ~ M l (42) 
Vi ~ M l , S ± ~ M 1 (43) 

for example, the first stemming from the Si = SiSqSi term for the scalar piece and from the SiSoVi, SiVqSi terms 
for the vector piece. Note that the system has no fixed-point solution with S\ ~ M n with n > 1 (the SiSqSi term 
then runs away upon iteration) nor with V\ ~ M n with n > 1. This still leaves us with the situation that, for a 
self-consistent solution both the scalar and vector pieces grow with the quark mass. As will be seen later on in section 
IVIII this yields an enhanced string-tension between heavy-quarks that grows quadratically with M. 5 

Let us now back-feed these solutions to the quark-propagator DSE to ascertain whether any change in the per- 
turbative 1/M behavior of the propagator is to be expected. Since the inverse propagator is the amputated Green's 
function with two quarks and no gluons, by analogy it is obvious that the 1/M counting for higher-point functions 
(with two quarks, n gluons) is the same. Consider first the vector part of the propagator. We have Vq > 0(1) because 
of the presence of the un-suppressed vector term of the inverse bare propagator /k on the right hand side of the 
DSE. Substituting now V\ = 1, M in the 7V0V1 structure of the loop we obtain M~ 2 or M~ 2 M = M _1 . Thus the 
vector part of the propagator remains bare in the heavy-quark limit. Consider now the scalar part of the propagator 
equation. This gives Sq > 0(M) from the bare mass term of the propagator. If we substitute S\ = M and Vq = 1 
into the self-energy, the relevant structures 7 So V\ and 7V0S1 are 0(M^ 1 ) and O(M ) respectively, and therefore 
sublcading compared to the bare mass term. Also, the solution with Vi = 1 yields a power M°. We are thus left with 
the purely perturbative counting in the presence of a heavy quark mass scale, that is, 

V ~ M°, S ~ M 1 , (44) 

consistent with our assumption above. 

For the quark-gluon vertex we are thus left with all the three possibilities (|4"Tj) . (|4"2"|) and (|4"3")l . It is then a purely 
dynamical question which of these three possibilities are realized. As a matter of fact there is even a fourth possibility 
which we discuss now. 



2. Modified counting: M (0) 9^ m 

There is yet another logical possibility that solves the coupled system of propagator and vertex equations, and this 
is a propagator that is not suppressed by a power of the heavy quark current mass in the infrared, that is, with some 
degree of saturation in M{p 2 ) = B/A, a mass- independent constant. Substituting this in the vertex equation, it is 
obvious that self-consistency is reached by a quark-gluon vertex that is also mass- independent, i.e. we have 

Vi - M° , Si - M°. (45) 

This case has interesting consequences. As will be discussed later on in section IViTl the string tension is then mass- 
independent. However, in the Schrodingcr equation one cannot take the infinite-quark mass limit and simultaneously 
study the two quarks at infinite distance. The infinitely heavy-quark mass limit only makes sense at intermediate 
momentum where the propagator is still affected (suppressed) by the heavy quark mass. 



We presently ignore how the lattice data would rule out such a scenario. It is far from intuitive since lattice practitioners often fix the 
scale to the string tension, now scaling with the quark mass, and in the heavy mass limit all other splittings they could measure in 
heavy quarkonium also scale with the mass. Happily experimental data, as seen later, seems to rule this solution out. 
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As mentioned before, the four cases (|4Tj) . (|42|) . (|43|) and (|45|) are all selfconsistent and it is a question of dynamics 
which one of them is realized. It may even be, that in different regions of the quark mass different scenarios are 
realized, i.e. one may e.g. encounter the scaling law (|45|) for physical masses up to the bottom quark mass, whereas 
the naive perturbative scaling (|41[) is realized for quark masses as large as the top quark. 

In fact, as we will see in section lPVl the scenario (|45|) is favored by our numerics. The vertex is only mildly dependent 
on the current quark mass for masses at least up to the bottom quark. 

However, this behavior does not persist for larger masses. Our numerical calculations show that for masses in the 
range of 20 — 40 GeV the scaling goes over into the perturbative one given by eq. (jlTj) . 



I. The Slavnov-Taylor identity 

In this subsection we check that the infrared power-counting that we have devised satisfies the pertinent Slavnov- 
Taylor identity for the quark and gluon vertex. Of course, since our analysis largely hinged on dimensional reasoning 
one does not expect to find inconsistencies here. 

In the following we assume all momenta to be small and proportional to each other, such that our infrared counting 
rules can be applied. 

The Slavnov-Taylor identity for the quark-gluon vertex is explicitly given and derived in (6~H . One needs to introduce 
a non-standard vertex H a (k,q) that corresponds to a ghost-to quark off-diagonal self-energy. This of course is an 
auxiliary device that violates all of fermion number, color and momentum conservation. One can think of it as a 
quark-ghost source. This vertex can be obtained in terms of an equivalent "bare vertex" gT a via an inhomogencous 
Bethe-Salpeter equation 

H a (k,q)=gT a + B a (k,q) (46) 
with the loop integral coupling it to the ghost-quark scattering kernel. 




-B a {k,q) . 



The momenta k and q correspond respectively to the incoming ghost and quark legs. 

The function H is infrared finite within our power counting. To sec this, expand the quark-ghost scattering kernel 
in terms of a skeleton graph and count the anomalous dimensions of the lowest order using the general counting rule 
([5]) and the infrared leading l/A/-terms in the quark-propagator. One arrives then at a total infrared anomalous 
dimension of k — 1. This power is precisely canceled by the — k — 1 from the ghost propagator in the diagram above 
and the +2 from the loop integration measure. Hence the total exponent is and the function H takes a finite 
constant value in the infrared. 

The Slavnov-Taylor identity for the quark-gluon vertex reads 

p£r°(pi,Pa) x G-\pl) = H a ( Pl ,p 2 )S- 1 (p 2 ) - S- 1 (p 1 )H a (p 1 ,p 2 ) (47) 

in terms of the inverse quark propagator S 1-1 and inverse ghost self-energy G~ 1 (note this is written down as 1 + fe(p§) 
in the original reference (64j but there it was assumed to be infrared finite as opposed to infrared suppressed). Since 
the quark propagator in the chirally symmetric broken phase is a constant and we just argued that H is also an 
infrared constant, so is the right hand side of the STL 

The left hand side power-law is easy to count: —1/2 — k from the quark-gluon vertex, 1/2 from the explicit gluon 
momentum, and +k from the inverse ghost dressing G _1 precisely cancel out to yield 0. Thus our infrared power- 
counting is consistent with the requisite gauge invariance. Note also that the argument is independent of whether 
the quark momentum is small or not, so it applies to either soft or uniform divergence counting, as long as the same 
counting is used for the quark-gluon vertex and the quark-ghost scattering kernel. 

Finally we discuss the case of the chirally symmetric solution. Here the H-function is still constant since to lowest 
order in the skeleton expansion there are two powers of (p 2 ) -1 ^ 2 less in the quark-gluon vertices but also two powers 
of 1/p more in the quark propagators. The inverse propagator S on the right hand side of the STI is now given 
by i fr, whereas the vertex dressing on the left hand side cancels the contribution of the ghost dressing function. So 
both the left and the right hand side of the equation are proportional to p and the matching still works. 

It is also easy to see how at least the perturbative and modified mass- countings for the quark and gluon vertex, 
for heavy-quark masses, are consistent with the Slavnov-Taylor identity. In the modified counting, neither propagator 
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nor quark-gluon vertex are dependent on the quark mass, therefore there is no running in cither LHS or RHS of the 
STL The perturbative counting can be understood from the Abelian WTI in QED 

i(fc„ - k'^T^k') = S-\k) - S- 1 ^') (48) 

where for fixed k, k' and M — > oo the RHS appears to have a vector part of the fermion propagator Vo oc M°, and a 
scalar part So oc M 1 ; however, neglecting momentum dependence, the scalar part on the RHS cancels (M — M = 0) 
and since the LHS stays mass-independent, because 7 M is mass-independent, the mass independence of both sides 
coincides. Coming now to the QCD STI for the qqg vertex in perturbation theory, in this counting the H a kernel 
becomes a mass- independent constant at large mass (only the gT a bare term survives) and the M — M cancellation 
happens again on the RHS. Thus at least the perturbative and modified large-mass countings are seen to be consistent 
with the STI. For other mass countings the situation is more complicated and we will not discuss them further. 

IV. NUMERICAL RESULTS FOR THE QUARK-GLUON VERTEX 

Having discussed the infrared structure of the propagator and vertex DSEs analytically, it is now time to underpin 
our findings with numerical studies. To this end we solve the coupled system of the DSE for the quark propagator 
and the one for the quark-gluon vertex in the infrared preserving approximation discussed around fig. [5J For the 
convenience of the reader we display the resulting system again diagrammatically in fig. 1111 In the following we first 
discuss our numerical techniques and then present our results. 

A. Renormalization and numerical methods 

In our numerical calculations we use solutions of the coupled ghost and gluon propagator DSEs as input, cf. appendix 
IC 31 These have been renormalized at the scale /i 2 = 170GeV 2 , which we also use to renormalize the quark-gluon 
vertex and the quark propagator. The corresponding renormalization conditions for the quark wave function Zf(p 2 ) 
and the vertex r„(p 2 ,p2,p 2 ) are 

Z f {f) = 1 

Ai(/! 2 ,2// 2 ,3/! 3 ) = 1. (49) 

Note that only Ai appears in this condition, since the accompanying tensor j v is the only divergent part of the 
quark-gluon vertex; all other dressing functions A2..4 and ri..g are UV-finite. Note also that we renormalize the vertex 
at the totally asymmetric momentum slice (cf. appendix IA 2[) . which is convenient for our numerical setup. The 
conditions (|49|) then unambiguously determine the renormalization constants Z2 and Z\f for the quark propagator 
and the quark-gluon vertex. 




FIG. 11: The DSE for the quark propagator and the quark-gluon vertex as used in our numerical investigation. The open 
circles in the diagram on the right hand side of the vertex-DSE indicate that we use fit functions for these internal vertices, 
whereas the quark-gluon vertex in the quark-DSE is the one calculated from the vertex-DSE. 
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The loop integral on the right hand side of the DSE is treated employing standard techniques. We perform all four 
integrals numerically in hypcr-sphcrical coordinates using Gauss-Legendre integration; no angular approximations are 
done. The radial integral is treated on a logarithmic integration grid with explicit infrared and ultraviolet momentum 
cutoffs. We verified that our results are independent of these cutoffs. 

The methods to project the right hand side of the vertex equation onto the various tensor structures of the vertex 
are described in some detail in appendix[G] We evaluate the vertex from the vertex-DSE on a three-dimensional grid, 
which is adapted to the quark DSE. To this end we use the propagator fits (|C9|) and (|C10[) as our start guess for the 
quark propagator. Upon backfeeding the vertex into the quark DSE, this equation is iterated until convergence. The 
resulting quark propagator is then again back-fecded into the vertex-DSE and so on until complete convergence of 
the coupled system of equations is achieved. 

Due to the considerable numerical complexity we have not yet found a way to completely back-feed the vertex into 
its own DSE. We therefore employ only parts of its tensor structure for the two internal vertices in the non-Abelian 
diagram. These are represented by fit-functions which are matched against the calculated vertex in the asymmetric 
kinematical momentum slice. From our infrared analysis in subsection ((III Dj) we saw that in principle it is sufficient 
to back-feed the most important vector and scalar components into the equation to generate selfconsistent infrared 
physics (at least on a qualitative basis). We therefore employ the two structures 

r M (pi,l>2) = JfiM(pi,P2) - i(pi +P2) A1 A 3 (pi,p 2 ) (50) 
for the internal quark-gluon vertices in the non-Abelian diagram. Here we will use the fit forms 

\i(pi,P2)=f a (x)gT(x) , A 3 (pi,p 2 )= 7= 1 = f n (x)g?(x) (51) 
with x = p\ + p\ + p\ . These include the proper asymptotics in the IR and UV 

9/44 / , \ «3 / j \ 2 

. »rw = ( . (52) 



X 



^2 + x ni J \d 3 + x / 

For the three-gluon vertex we use the tree-level tensor structure and represent the corresponding dressing function by 

-3k / j r x -I \ 17/44 

+ d 2 log — + 1 
cti + x J \ d\ + x «i 



r°' 3 ( Pl ,P2) = ( ) ( d,-fi- + d 2 \o g 



(53) 



which reflects the strong (p 2 ) _3K singularity known from Yang-Mills theory, eqs.([5|). (fTU|) and a logarithmic ultraviolet 
tail with coefficients from resummed perturbation theory. In all fits the momentum argument is given by the sum of 
the three squared momenta from the external legs. The fit forms thus represent overall divergences only and do not 
take care of possible soft gluon divergencies. We will modify the fits accordingly when we investigate this issue in 
subsection IIVDI 

In general, backfeeding only parts of the vertex into the vertex equation is somewhat problematic. Since most of the 
tensor structures of the vertex contribute considerably to the loop integral of the non-Abelian diagram, omitting some 
of them leads to reduced interaction strength in the integral. As can be seen from our power counting in subsection 
IIIIDI this fact does not affect the qualitative behavior of the resulting vertex dressing functions (i.e. their infrared 
exponents), but does influence the coefficients of the resulting power laws. As a result one faces a dilemma: 

(i) One can insist that the input fit functions match both the output infrared exponents AND their coefficients. 
We found that this is only possible if our input for the three-gluon vertex is v ery strong, i.e. much stronger 
than one would expect from corresponding DSEs for the three-gluon vertex [54| . The missing interaction 
strength, presumably from the neglected ggg vertex pieces is then simulated by the strong three-gluon interaction. 
This setup allows one to systematically compare the resulting quark-gluon vertices for different quark masses. 
However, as we will see, it leads to unnaturally large scales, which cannot be true. 

(ii) If we instead insist on correct physical scales, measured e.g. by the chiral condensate from the quark propagator 
in the chiral limit, the missing interaction strength in the non-Abelian diagram does not permit a match to the 
scales from the input fits with the resulting vector and scalar vertex pieces. The missing interaction strength 
from the neglected vertex pieces in the loop diagram is then taken over by the input dressing functions, which 
are stronger than the calculated ones. 

Naturally this quantitative problem would not exist if we were able to back-feed the complete vertex into the vertex- 
DSE. Postponing this considerable numerical effort for future work, we will use a combination of the possibilities (i) 
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and (ii) here. We first investigate the chiral limit using the procedure (ii) and thus producing phcnomcnologically 
acceptable scales in the resulting vertex and quark propagator. For our quark mass study, however, we have to resort 
to the procedure (i). Since we are dealing with qualitative statements only, this compromise seems acceptable to us 
and we believe that the results we obtain are meaningful. 

B. The chiral limit 

Let us first have a look at the chiral limit and the corresponding structure of the quark-gluon vertex. The input 
parameters in our calculation for the fits f51"j) and |53|) are given in table [ill Here the ultraviolet behavior of the three- 
gluon vertex is controlled by the parameter cfo, which is chosen such that the vertex is normalized to the tree- level 
expression at the renormalization point /i 2 = 170GeV 2 . The corresponding ultraviolet coefficients d 2 , d^, n\ and n 2 
of the input functions Ai and A3 are chosen such that the calculated vertex dressings are matched by the input. The 
remaining infrared coefficients d\ in Ai,3 are chosen such that the relative strength of the corresponding calculated 
dressings is matched. Finally, the parameters d\ and d% in T 0,3 are chosen such that the chiral condensate of the 
quark propagator reproduces the phcnomcnological value. 
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TABLE II: Input parameters for the fit functions (|51|l and (|53|) in the chiral limit. 

The resulting input and output vertex dressings for Ai and A3 are shown in fig. [12] together with the resulting quark 
propagator and further vertex dressing functions. 

In the top right corner we compare the internal representations of the vector (7^) and the scalar (p^) piece of the 
quark-gluon vertex with the calculated components. The most important observation here is the perfect matching 
of the generated infrared exponent (p 2 ) -1 / 2 ~ K of both the Ai and A3 components with the infrared exponents from 
the input fits. In fact the very same power law is developed by all twelve tensor structures A1..4 and Ti..g- One then 
expects that backfeeding additional structures besides Ai and A3 would just reproduce the same infrared exponents 
as seen from fig. [T3]with backfeeding Ai and A3 only. We have explicitly checked that this is indeed the case. Thus 
the numerical result confirms that all dimensionless vertex dressings share the same infrared divergence 5 . 

For the quark mass function M(p) we observe the characteristic rise towards small momenta that is typical for 
dynamical chiral symmetry breaking. Quite unexpectedly we also observe a turnover at very small momenta, which 
to our knowledge has not been seen yet in any other truncation scheme. Whether this turnover is a genuine effect 
independent of the truncation of the vertex-DSE remains to be investigated. Also for the quark wave function Z_r(p) 
we observe a stronger decrease in the infrared than has been observed in previous truncation schemes (see e.g. [9(). 

As concerns the scale we did not match the input and output scales here but pursued scenario (ii) and seeked to 
generate the quark-gluon vertex at a 'physical' scale such that the chiral condensate has the correct magnitude. This 
condensate can be extracted from the quark propagator in the chiral limit according to 

- (**>(p) := Z 2 Z m N c tr J A^Soiq 2 ), (54) 

where the trace is over Dirac indices, So is the quark propagator in the chiral limit and Z m is the quark mass 
renormalization factor that can also be determined from the quark-DSE. We obtain 

(**)^' =0 (2 GeV) = (-260) 3 McV, (55) 



Ai 


A 2 


A3 


A 4 


n 


T2 




TA 


TT, 


T6 


T7 


T8 


1 


Ap( 


2pi 


2pi 


A 


2p\ 


P i 


Apt 


Pi 


2pi 


2 P { 


2p( 



TABLE III: Dimensionful factors with which the vertex dressings are multiplied to obtain dimensionless quantities. 



In the notation of appendix [D] and ref. |57ll this amounts to f3\ = Pi = Ps = Pi = 



-1/2 - K. 
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FIG. 12: Top left corner: the chiral quark propagator on a linear momentum scale. Note that the maximum seen in the 
quark mass function is not a strict prediction of our truncation and may or may not be a truncation artefact. Top right corner: 
internal representations of the vector and scalar vertex dressings (dashed lines) compared to the corresponding calculated vertex 
pieces (full lines) . Bottom left corner: The six large tensor components of the vertex. Bottom right corner: the six small tensor 
components of the vertex. Note the different momentum-scale as compared to the plot of the large components. All vertex 
dressings have been multiplied by appropriate powers of momenta, given in table II, such that they correspond to a normalized 
version of the basis eq. (2) and can be compared quantitatively. The dressings are plotted versus p = pi and and evaluated in 
the asymmetric momentum slice with p\ — ?>p\ and p\ = 2p\. 



where we have used A-p-g = 225 MeV. This value is in agreement with the values reported in the literature. 

In the bottom line of fig. [T^] we compare the strengths of the different structures of the vertex. These are multiplied 
by appropriate factors to make them dimensionlcss; the factors are given in tabic UlTl Clearly sonic structures arc more 
dominant in terms of infrared strength (i.e. coefficients of the universal power law (p 2 )~ 1 ^ 2 ~ K ), the most important 
pieces being Ai, A3 and T5 7 . There is no uniform behavior of the various vertex dressings but some are strictly 
positive and monotonic, some are strictly negative and some change sign in the infrared. This behavior may or may 
not be an artifact of the chosen truncation scheme of the quark-gluon vertex DSE and will be investigated further in 
future work. Of all structures perhaps the most interesting one is the A3 piece, since it represents the leading chiral 
symmetry breaking component of the vertex. Its considerable infrared strength signals a large amount of dynamical 
chiral symmetry breaking not only in the quark propagator but also in the vertex structures. This is consistent with 
and in fact corroborates the strategy and results of our infrared analysis in subsection IIIIDI and appendix |Dj 



7 This has already been noted in fl5| , |16| . 
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FIG. 13: The vertex dressing function Ai(pi,p2) at fixed angle arccos(pi ■ f>2) = 0.1. 



A two-dimensional plot of the dressing function Ai(pi,p2) at fixed angle arccos(pi ■ P2) = 0.1 is shown on in fig. 1131 
One clearly observes, that the infrared divergence of Ai occurs almost uniformly in the p\ — p|-planc. This is in 
qualitative agreement with lattice calculations of the vertex [TlT ] at those momenta available on the lattice so far. 
Unfortunately the volumes of these lattice simulations are by far not large enough to be able to confirm or refute the 
divergence of the vertex in the deep infrared. 

Next we determine the running coupling from the quark-gluon vertex. For the convenience of the reader we repeat 
its definition here, (c.f. subsection IIII G[) : 

a qa {p 2 ) = ^-\ 2 (p 2 )Z 2 {p 2 )Z{p 2 ). (56) 

Although the right hand side of this definition is renormalization point invariant, it still depends on the renormalization 
scheme. The scheme we are using, eq. (|49]l . is adapted to the asymmetric kinematical momentum slice. We call 
this an 'asymmetric MOM-schcme'. The resulting running coupling as numerically calculated in this scheme is 
displayed in fig. 1141 Clearly one sees the usual logarithmic running of the coupling in the ultraviolet, with a value 
a MOM-asyrnm( M 2^ _ q 118 at tlle mags Q f the z_b 0S on. (Note that the exact agreement of this result with the 

experimentally determined coupling in the MS — scheme is incidental.) Towards lower momenta we see a gradual 
increase of the coupling, where the nonperturbative region a MOM-asymm y. \ \ s reached at about p 2 = 0.75GeV 2 . 
The coupling then continually grows, seems to reach a plateau (corresponding to a similar plateau in Ai) and then 
enters the scaling region in the infrared. In this region, below p 2 = (lMcV) 2 we see a 1/p 2 singularity in agreement 
with our analytical analysis in subsection IIII Gl 

Finally we investigate the possibility of positivity violations in the quark propagator. These are related to the 
analytic properties of the quark propagator which can in part be read off from the corresponding Schwinger function 

f • f d A p 

<r(t) = J d 3 x J -j—-^ exp(ip ■ x)a S y{p 2 ) 1 (57) 

where os.v are the scalar and the vector parts, respectively, of the dressed quark propagator. (This method has a long 
history, see e^q. [HD, [66|, [6?], [68| and references therein) . According to the Osterwalder-Schrader axioms of Euclidean 
field theory [69| , this function has to be positive to allow for asymptotic states in the physical sector of the state space 
of QCD. Conversely, positivity violations in the Schwinger function show that the corresponding asymptotic states (if 
present) belong to the unphysical part of the state space. Thus positivity violations constitute a sufficient condition 
for confinement. Our results for the Schwinger function of the chiral limit quark propagator are shown in the right 
diagram of fig. [14] Shown is the absolute value of the Schwinger function on a logarithmic scale. The thick black 
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FIG. 14: Left diagram: The running coupling from the quark-gluon vertex in the asymmetric MOM-scheme with /i 2 = p\ = 
P2/2 = p|/3. Infrared slavery in the form of a power-law is seen at extremely low-momenta, whereas the coupling is quite flat 
at intermediate to low-momenta, just as needed for conformal arguments (65[ . Right diagram: absolute value of the Schwinger 
function of the scalar part of the quark propagator in the chiral limit. Clearly visible are positivity violations at a scale of 
about 5/GeV ~ lfm. 



curve denotes our numerical result for the Fourier transform (|57j) . Clearly visible is a cusp in the Schwinger function 
at about iRi5 GcV -1 f« 1 fm. Certainly, this is a typical scale where quark confinement is expected to occur. 

A possible analytical structure for such a positivity violating quark propagator has been suggested in [70| . It reads 

ait) = |6oexp(-M)cos(6 2 i + & 3 )| > ( 58 ) 

which corresponds to a pair of complex conjugate poles of the propagator in the timelike momentum plane. These 
poles correspond to a 'quark mass' given by m = b\ ± 162- In our case we can fit this form either to the small time 
behavior of the numerical result or to the region between the first two cusps, as shown in the diagram. The results 
are then m = 305(5) ± i 453(10) MeV for 'Fit A' and m = 351(5) ± i 374(10) MeV for 'Fit B'. Here the errors reflect 
variations in the bounds of the fitted regions. 8 



C. Mass dependence of the coupled system 

We now investigate the dependence of the quark-gluon vertex on the current quark mass. To compare the results 
for different quark masses in a meaningful way we match the input and output scales of the vector Ai and scalar A3 
components of the vertex in a given momentum configuration; since the results are qualitatively similar for different 
kinematical slices, they are most easily presented using the asymmetric momentum slice (c.f. subsection IA2| . This 
choice is not unique, since we saw in a previous subsection that the vertex is not uniform with the angle pi ■ p2 in 
the infrared, whereas our input fits are. However, we checked that different choices would not affect the qualitative 
statements of this subsection, so the asymmetric slice can be chosen without loss of generality. 

The input parameters in our calculation for the fits (|5ip and (|53[) and the current quark mass are given in table [TV] 
Again the ultraviolet scale of the three-gluon vertex is matched such that the vertex goes bare at the renormalization 
point fj 2 = 170GeV 2 . The infrared strength of this vertex is now considerably stronger than in the previous sections: 



This result demands a comment. In 15311 two of us showed that the quark propagator is positive definite if an ansatz for the quark-gluon 
vertex is used which contains an Abelian part, satisfying the Abclian Ward-Takahashi identity, and an overall non-Abelian dressing 
function multiplying all tensors uniformly. The resulting DSE for the quark propagator looks similar in structure to the corresponding 
one for the fermion propagator in QED. Consequently one obtains an analytical structure of a singularity at the fermion mass with a 
branch cut attached, i.e. exactly the structure expected for a physical electron. It was then left open in ref. [53ll whether the exact 
quark-gluon interaction of QCD has such a structure, leading to a positive definite quark confined by another mechanism, or whether the 
quark-gluon interaction is sufficiently different from this ansatz so as to allow for positivity violations. Our results here point in exactly 
this direction. However, since the relative strength of the vertex components are somewhat truncation dependent further research is 
necessary to confirm this result. 
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TABLE IV: Input parameters for the fit functions (|51[) and (|53[) and the current quark mass. The parameters are chosen such 
that the input fits to Ai and A3 match the calculated dressing as closely as possible in the asymmetric momentum slice. 



u/d-quark 




10 
P [GeV] 



10 



10" 



10"" 



b-quark 




— \ 


-fit 


— x 3 


-fit 


— \ 




— h 









lJ I I I I 



10" 



10" 

P [GeV] 



10 



10 



10 



FIG. 15: Internal representations of the vector and scalar vertex dressings (dashed lines) compared to the corresponding 
calculated vertex pieces (full lines). Shown is the asymmetric momentum slice with p% = Zp\ and p\ = 2p\. 



the vertex now effectively takes over the missing strength from the vertex structures which are not back-feed into the 
loop integral. The remaining parameters in the Ai,3 input fits are then chosen to match their calculated counterparts. 
Since we are matching in a particular kinematic slice only, it does not make sense to seek for perfect matching here; 
we therefore concentrated on the infrared power law coefficients and the ultraviolet tails and left some discrepancies 
at the intermediate momenta. These discrepancies are irrelevant for the points we wish to make here. 

The resulting input and output vertex dressings for the two extreme cases (u/d and bottom quark) are shown in 
fig. [121 One observes the following: 

• In the infrared the calculated vertex dressings reproduce the power laws from the internal fits regardless of the 
quark mass. This observation nicely confirms our infrared analysis from subsection IIIIDI and underlines the 
presence of universal power laws with mass independent exponents in the infrared. 

• We now come back to a point already discussed at the end of section IIII Bl The presence of these mass 
independent infrared anomalous dimensions also shows that the potentially small scale Aqcd/M is not small 
enough (at least for physical quark masses) to disturb the power counting in the infrared. In fact we do not see 
any effects from this scale in our numerical results at all. This again justifies a posteriori our infrared analysis 
in sections IIII Dl and [D] 

• Comparing the coefficients of these power laws, given in table |V] we see a much smaller dependence on the quark 
mass than could have been anticipated from the naive analysis in subsection IIII Hi Assuming some systematic 
errors one could interpret the results for both, the vector Ai and the scalar A3 dressing function as constant 
if not slightly decreasing with the quark mass. The value of the quark mass function at zero momentum is 
only slightly increasing with the current quark mass. Clearly, this pattern does not match any of the three 
possibilities found in our analytical analysis with 'naive' counting in subsection IIII H l\ Where does the naive 
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FIG. 16: (Color online) The fully converged quark propagator from the coupled system of vertex and propagator DSEs for 
the different current quark masses (from bottom to top) m u / d (p?) = 3MeV, m s (fi 2 ) = 50MeV, m c (fi 2 ) = 1200 MeV and 
rrib(fi 2 ) = 4200 MeV at the renormalization point p 2 = 170 GeV. 



counting of quark masses go wrong? The basic assumption of the naive counting is an increase of the infrared 
mass M = M(p 2 -C Aq CD ) roughly proportional to the current quark mass. This is, however, not the case in 
our solutions as can be seen from fig. [TSJ This could be in agreement with subsection IIII H 21 Whereas the 
ultraviolet behavior of the quark mass functions indeed is proportional to the current quark mass (as strictly 
demanded from perturbation theory), the quark mass in the infrared is not. Whereas the current quark mass 
varies by a factor of O(10 3 ) from the up- to the bottom-quark mass, the infrared mass M varies only by roughly 
a factor of two. Together with the variation of Zf{p 2 ) in the infrared this mass dependence explains the result 
for the mass (in-)dependence of the quark-gluon vertex observed above. 

• We also did calculations for very heavy quark masses in the range of 20 — 40 GeV. In this region, and presumably 
also for heavier quarks, the naive scaling again sets in, i.e. M (0) cx m and we see perturbative counting behavior 
for the vertex. It therefore seems that the region of physical quark masses is clearly distinct. 

These two last observations are somewhat unexpected, since they deviate from the usual behavior found in the 
literature. It remains to be seen whether this effect persists if the quark-gluon vertex is back-coupled self-consistcntly 
into its own DSE. Nevertheless it is interesting to speculate of the possible consequences of this result. This will be 
done in section IVIII where we derive some qualitative aspects of the string tension related to an infrared divergent 
quark-gluon vertex. 



D. Soft-gluon divergence 



In the last subsection we showed that indeed the quark-gluon vertex as calculated from our vertcx-DSE has a strong 
divergence like (p 2 ) _1 / 2_K in the infrared once all external scales go to zero, i.e. p 2 ~ p\ ~ p\ ~ p\ — > 0. However, in 
our infrared analysis, subsection IIII El we also discovered an additional selfconsistent infrared singularity when only 
the external gluon momentum goes to zero, i.e p 2 ~ p 2 — ► 0, but p\ and p\ are kept constant. We now confirm this 
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Ai[GeV 1/2+K ] 


3.95 


3.60 
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2.73 


X 3 [GeV 1/2+K ] 


8.70 


7.97 


4.44 


6.08 



TABLE V: The infrared coefficients for the dressing functions Ai,3 for different current quark masses m(n 2 ) at the renormal- 
ization point /i 2 = 170 GeV. Also given is the (approximate) value of the mass function at zero momentum. 
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FIG. 17: Input and output of a coupled quark-propagator (totally self-consistent) to quark-gluon vertex (one loop iteration 
with internal dressings) showing that the soft-singularity is indeed a property of the self-consistent vertex equation with the 
non-Abelian diagram. 



analytical result also numerically. To this end we change the internal dressing functions (|5ip to the following forms: 
Ai(pi,j*) = / IR (p!)<?r(z) , A 3 (pi, P2 )= n 1 = m)g^(x) (59) 

which have additional strength in the soft gluon point, i.e. they go like ((P3) 2 ) ^ 2 where P3 = P2 — Pi is the gluon 
momentum. Similar as in eqs. (j5T]) we have x = p\ + p\ + p\ for the ultraviolet part of our fits. The parameters 
are unchanged, i.e. given in table [II] For the quark propagator we use the chiral result of the last subsection. In 
fig. [17] we plot a calculation based on these input forms. Shown is the output dressing function Ai in the kinematical 
section with p\ ■ p% = 0. We fix p\ = 0.24 GeV and vary p\. The soft gluon point p\ = is then reached when 
p\=p\. As can be seen from the plot we clearly encounter a singularity in this limit. From a corresponding log- log 

plot we extracted the exponent of the power law and found ((pi — P2) 2 ) 12 ^ to good accuracy in agreement with 
our infrared analysis. 

This observation has several important consequences. First, the limit p\ —> is important for the confining potential, 
as detailed later in section IVIII Second, the soft-gluon endpoint p§ = is a preferred point chosen in several lattice 
calculations [ll|, G3- However, we just saw that the vertex is singular at precisely this kinematical section. This then 
begs the question: how can lattice calculations obtain results at this point. 

If it is true that a soft-gluon power law develops, then at fixed gluon momentum pz = one should obtain an 
infinite result no matter what quark momentum is chosen. Thus, a computation in the soft-gluon kinematic section is 
bound to be infinite. On the lattice, however, this infinity is regulated and one obtains a finite answer for all momenta 
P 2 — Pi = p\ a t p\ = 0. This finite answer should then vary with the volume and lattice spacing, provided there 
is indeed a soft-gluon singularity in the infinite- volume/continuum limit. We therefore propose to study the scaling 
properties of the vertex in the soft-gluon kinematical section. There is also one further subtile point here: all scaling 
effects may be wiped out if the vertex is also renormalized in this section (as partly done in [ll|, (HI). Therefore it is 
recommended to change the renormalization point to a non-divergent section, as for example our totally asymmetric 
kinematic section p\ = l/2p?, = l/3p§. Then the graph of Ai in the soft-gluon kinematic section is not forced to pass 
by one, and the sensitivity to varying the grid is exposed. 

This exercise is carried out in fig. 1181 where we indeed find a large sensitivity of the results to the grid in the 
presence of our soft-gluon divergence. 
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Detecting a soft-gluon divergence in the qqg vertex 

(Plot symmetric slice, renormalize at asymmetric slice) 
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FIG. 18: Black solid line: input Ai dressing for the internal vertex with a soft-gluon singularity, plotted here against the 
momentum of the gluon in a safe kinematic section. The calculation also employs a A3 dressing, not shown. Red dashed line: 
the output Ai dressing after one iteration of the non-Abelian loop with the Curtis-Pennington quark propagator as input. The 
grid has 60 points and is cut-off in the infrared at about 1 MeV (therefore the loop integral too). The plot is in the soft-gluon 
kinematic section where the analytical result is divergent. Therefore, the curve is an artifact of the grid regularization. Blue 
dotted line: same but with an 80-point grid and a lower cutoff of 0.01 MeV. As seen, the grid sensitivity is large when the 
renormalization point is not imposed on this kinematic section, forcing the curves to pass by 1. This spurious effect might be 
at work in existing lattice calculations. 



E. Chiral symmetric solution 



Finally we present numerical solutions for the case of unbroken chiral symmetry. This means we solve the coupled 
system of quark propagator DSE and vertex-DSE employing chirally symmetric tensor structures only. As detailed 
in subsection IIII Fl the quark propagator is then given by 



S X (P) 



if> 
V 1 



Zf(p 2 )- 



(60) 



In general kinematics the corresponding chirally symmetric quark-gluon vertex has the tensor structures L\, L2, 
T2, T3, Tq and T% given in eqs. (|A2[) with dressing functions A1.2 and T2.3 i 6.8- All other dressing functions vanish 
identically. This is reproduced by our numerical solutions. 

Since A3 is among the vanishing dressing functions we only work with Ai in the internal quark-gluon vertices. The 
corresponding fit function is given by 



Ai(pi,_p 2 ) 



di 



d 2 log 



x 



-9/44 



(61) 



Note that the infrared exponent is changed from —1/2 — n, as in eq. I|5ip . to — k and we have introduced an additional 
parameter d^. The input for the three-gluon vertex is unchanged compared to (|53|) . Again we matched the input 
and output values for Ai in the asymmetric momentum slice p 2 — pf = = p§/3. The resulting values for the 
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Ai 


0.10 


0.20 


0.02 


■p0,3 


1.00 


0.40 


1.00 



TABLE VI: Input parameters for the fit functions (|5ip and (|53p in the chiral limit. 



29 




FIG. 19: All diagrams show results for the quark-gluon vertex, the quark propagator and the running coupling in the case 
where chiral symmetry is forced to be unbroken. Top left corner: Internal and external (calculated) dressing function Ai in the 
asymmetric momentum slice with p§ = 3pf and p\ = 2p\ . Top right corner: The quark propagator function Zj(p 2 ). Bottom left 
corner: The six chirally symmetric tensor components of the quark-gluon vertex. Bottom right corner: The running coupling 
for the chirally symmetric case. Note that within numerical errors the coupling approaches a fixed point in the infrared in 
agreement with our analytical analysis. The bump at mid-momentum is model dependent and potentially an artifact due to 
scale mismatch between the gluon propagator and the vertex functions. 



parameters are given in table IIVEI Here the parameters for the three-gluon vertex allowed for the 'natural' values 
d\ =(#3 = 1. These set the scale of the results. Certainly this scale is arbitrary, since we have no physical situation 
to compare with. 

Our results for the vertex and the quark propagator are displayed in fig. 1191 In the top left corner we show the 
matching between the internal fit form for Ai and our calculated result for the asymmetric momentum slice. The 
overall matching is very good, the slight deviations in the mid-momentum regime do not affect the point we would like 
to make: in the infrared there is a clear matching between the internal and external power laws once again justifying 
our analytical analysis from subscction llll Fl The vertex indeed diverges in this case with (p 2 )~ K . The resulting quark 
dressing function Zf(p 2 ), displayed in the top right corner, shows a somewhat unexpected behavior. It starts with the 
expected perturbative running in the ultraviolet momentum region. Then it bends down towards smaller momenta in 
the mid-momentum region similar to the chirally broken case. It then, however, starts rising again and settles for a 
finite infrared value larger than one. Taken at face value such a behavior forbids a positive definite spectral function 
[7lj . i.e. such a quark propagator does not describe a physical quark. Since this behavior may or may not be an 
artifact of our truncation it creates no problem per se, since the situation we are investigating in this subsection has 
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FIG. 20: The four-quark 1PI Green's function and the first terms of its skeleton expansion 



no counterpart in the real world. For our purposes it is sufficient to note that the quark dressing function Zf(p 2 ) 
indeed goes to a constant in the infrared, thus confirming our central assumption of the infrared analysis in subsection 



In the bottom left corner of fig. [TO] we show the dressing of the quark-gluon vertex as a function of the momentum. 
Similar to the chirally broken case we find Ai, T3 and ts among the large vertex components, whereas A2, T2 and tq 
are somewhat smaller. The running coupling in the chirally symmetric case is shown in the bottom right diagram 
of fig. 1191 In agreement with our infrared analysis we find something like a fixed point in the infrared. The scaling 
behavior is not as convincing as in the chirally broken case. To a large extent this may be to numerical uncertainties, 
which are much larger than in the chirally broken setup. We are working to improve this situation in the future. Note 
that the precise value of this fixed point depends on the input strength of the thrcc-gluon vertex and is therefore not 
a prediction of our calculation. It is the mere existence of this fixed point that will be important for our analysis of 
the quark-antiquark potential later on. 



In the previous two sections we found an infrared divergent, selfconsistent solution for the quark gluon vertex 
using both analytical and numerical methods. We now wish to investigate the consequences of such a behavior for 
the quark-quark scattering kernel, i.e. the Green's function that should contain the confining quark potential in 
the heavy quark limit of quenched QCD. We will first briefly discuss the relevance of chiral Ward-identities, then 
discuss the properties of four-quark Green's function in a skeleton expansion (cf. ref. [H3|) and finally show that the 
results derived from this expansion can also be directly justified from the Dyson-Schwingcr equation of the four-quark 
function. 



In the heavy quark limit chiral symmetry breaking is dominated by the effects of the large bare quark mass. Thus 
chiral Ward Identities, associated with effects from the dynamical breaking of chiral symmetry, play only a minor role. 
Certainly, the opposite is true when quarks with small masses are considered. Then, the axial vector Ward-Takahashi 
identity enforces an intimate relation between the scattering kernel and the gap equation. 

A systematic prescription to construct a scattering kernel that is consistent with a given quark gluon interaction 
has been given by Munczek [7^ |. Beyond-rainbow-ladder approaches have been demonstrated, for example in [73| . 
with semipcrturbative-type vertex constructions. The basic idea is to take functional derivatives of the quark self- 
energy with respect to the quark propagator. In appendix [Hi we detail a truncation scheme of the coupled system of 
Dyson-Schwingcr equations for the quark-gluon vertex, the quark-gluon scattering kernel and the quark propagator 
DSE that is amenable to Munczek's prescription. We then give a formal expression for a resulting Bethe-Salpeter 
kernel that is both chirally symmetric and confining. The construction is formally straightforward but leads to a 
complicated system of equations whose explicit treatment is far beyond the scope of this work. Instead we proceed 
discussing the heavy quark limit of the scattering kernel, where we expect a skeleton expansion to deliver meaningful 
results. 



We investigate the infrared behavior of the full quark-quark four-point function in the heavy quark limit and in 
the presence of only one dynamical external scale (the distance between a quark and an antiquark). This function 
together with the first two terms of a skeleton expansion are given in fig. 1201 As in subsection [D] we examine the 
infrared behavior of these graphs and perform a scaling analysis in terms of infrared anomalous dimensions. 



null 



V. THE QUARK-QUARK SCATTERING KERNEL 



A. Quark-quark scattering kernel and chiral Ward identities 



B. Quark-quark scattering kernel in a skeleton expansion 
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The first term of fig. [20] is nonperturbative one-gluon exchange, i.e. one-gluon exchange with a dressed propagator 
and two dressed quark-gluon vertices. Since we are working in the heavy quark limit with only one small external 
scale, i.e. the gluon momentum, we are probing the soft singularity of the quark gluon vertex (fc 2 ) _1 / 2_K . We have 
seen in subsection IIII El that this singularity appears in all tensor structures of the vertex. Thus we do not need to 
distinguish between vector and scalar components here. 

Together with the gluon propagator scaling (k 2 ) 2K ~ 1 we thus arrive at a total scaling 

terml - (fc 2 )" 2 (62) 

of the first term in the skeleton expansion of the four-quark kernel. Now observe how adding the second rung modifies 
the first term. For the quark propagator we have S — l ^ 2 f + ■jf —> ■jf in the infrared which leads to 



( fc 2)2+2«-l+2(-l/2-«) | ^ Z 2 
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term2 ~ terml x \ Z) ^ ^ + •■■ J ~ ( 63 ) 

for the second term. Here the dots represent sublcading contributions. The first factor of (k 2 ) 2 comes from the 
four-momentum integral, then we have (fc 2 ) 2K_1 from the extra dressed gluon propagator and two scaling factors 
(fc 2 )~ 1 /2-« f rom t h c 

quark-gluon vertices. We thus obtain the same scaling behavior of the second term than for 
the first term. With the methods of section IIIII it can be shown that this behavior persists for all other terms of the 
skeleton expansion. 

As a result we obtain for the four-quark kernel the power law H(k 2 ) ~ (k 2 )~ 2 . This is in perfect agreement with 
the 1/fc 4 behavior in the old infrared-slavery scenario. The crucial difference, however, is that in the former scenario 
the vertices have been assumed constant in the infrared, whereas the gluon was believed to be singular as 1/fc 4 . Here 
we have an infrared vanishing gluon compensated for by thc strength of two infrared divergent quark-gluon vertices. 

Then, the well-known relation 

y(r)= (2^/ d3fc ^ (fc ° = ' k)eikr;= (2^/ d3fc ? : ~ |r| (64) 

between the static four-quark function H(k = 0,k) and the quark potential V(r) therefore gives a linearly rising 
potential by naive dimensional arguments. We therefore find confinement in the chirally broken phase of quenched 
QCD. A more refined treatment, as described in leads to the same result. By construction already the first 
term in the skeleton expansion, i.e. nonperturbative one-gluon exchange displayed in fig. 1201 generates this result. 
Since the following terms in the expansion are equally enhanced in the infrared, the string tension will be built up by 
summing over an infinite number of diagrams. We will come back to this point below. 

It is instructive to also have a look at the chirally symmetric case. The corresponding infrared solutions for the quark 
propagator and quark-gluon vertex have been derived in subsection llll Fl and confirmed numerically in subscction llV El 
They correspond to a vertex scaling of (fc 2 )~ K and a gluon propagator which still goes like (/c 2 ) 2k ~ 1 . From the first 
term of the skeleton expansion for the quark-quark scattering kernel one obtains H(k 2 ) ~ 1/fc 2 . If one would derive 
a potential for such an interaction one finds 

nr)=pi F /^° = 0, k)e - = ^/^^ ~ (65) 

i.e. the well-known Coulomb potential. (Note, however, that the notion of a potential in chirally symmetric QCD 
is certainly not appropriate.) As we have seen in subsection IIII Gl also the resulting running coupling from the 
quark-gluon vertex is no longer diverging but goes to a fixed point in the infrared similar to the couplings from the 
Yang-Mills vertices. In a sense, the restoration of chiral symmetry is therefore directly linked with the disappearance 
of infrared slavery. This is one of the main results of this work. 



C. The four-quark Dyson-Schwinger equation 

We now proceed by justifying the above results directly from the Dyson-Schwinger equation for the four-quark 1PI 
Green's function, which is shown in fig. 1211 This equation describes the full amplitude that includes possible bound 
states, scattering as well as off-shell effects. In contrast to the skeleton expansion fig. [201 the diagrams on the right 
hand side of fig. [21] all involve a bare quark-gluon vertex. Since this vertex is infrared divergent it is therefore not 
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FIG. 21: The DSE for the four-quark 1PI Green's function. 




FIG. 22: Contribution to the skeleton expansion of the last term in the four-quark DSE fig. 1211 Soft and hard momenta are 
denoted by an s respectively h. 

obvious how both equations can give a similar infrared behavior for the four-quark function. Indeed, when considering 
only the graphs in the first line of fig.[2T]an infrared power counting similar to that in section Hill yields scaling laws 

(p 2 ) 2 ((p 2 )- K - 1/2 ) 3 (( P 2 ) 2K - 1 ) 2 =(/)- 3/2+K , (66) 

independently of the quark kinematics. Here, the second term arises from the three dressed quark-gluon vertices and 
the third one from the two gluon propagators, whereas the quark propagators proved to be scale invariant. These 
diagrams therefore cannot be responsible to generate a confining quark-antiquark potential. 

However, these diagrams are not the leading ones in the infrared. To see this, consider the last graph in the second 
line of the DSE, fig. [2TJ A skeleton expansion of this graph involves, among others, the two-loop graphs shown in 
fig. [5U In complete analogy to the analysis of the second ('t-channel') version of the DSE discussed in appendix [B] 
this two-loop integral has a kinematic region, shown on the left side (a) of fig. 1221 where the vertex correction given 
by the 'small' loop is dominated by hard loop momenta of the order of the quark mass and not by the infrared regime. 
This graph (a) represents the proper kinematics to study the confining properties of bound states. Here the momenta 
of the massive quarks in a meson are finite and do not scale to zero. However, the exchanged gluon momentum 
would become soft when the two quarks are sufficiently spatially separated. In this infrared limit the counting for the 
whole graph involves four dressed quark-gluon vertices in the soft-gluon limit p — ► where these vertices are divergent 
according to our analysis in section fill El This then yields 




i.e. a strong singularity in agreement with the one found in the skeleton expansion. As discussed above, this singularity 
has the potential to permanently confine quarks. 

Now consider the right graph (b) were the quark momenta are also soft and scale to zero. In Minkowski space 
this would correspond to the case of light quarks close to the on-shell limit, i.e. in some sense the opposite case to 
that where the quenched approximation we consider here is justified. Nevertheless let us study the power counting 
in this case. Naively the right graph shown in fig. [2U seems to yield an equally strong divergence than the preceding 
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one. However, there is no large quark mass here, so that the integral is dominated by small modes. Consequently 
all dressing functions are probed in their scaling region and it is the uniform limit of the quark-gluon vertices that is 
relevant here. The resulting power law is similar to the one of eq. (|66j) . i.e. not sufficiently divergent to trigger a linear 
rising potential. Comparing to the analysis of diagram (a) above, we see that (i) it is not the uniform infrared limit 
of the quark-gluon vertex that is responsible for quark confinement but the soft-gluon limit and (ii) it is impossible 
to confine light quarks close to the on-shell limit with this mechanism, whereas heavy quarks are confined. 

Unfortunately the four-quark DSE requires two-loop graphs and is far too involved to be useful for any explicit 
analysis of heavy bound states. Therefore in practical calculations one has to resort to the skeleton expansion given 
in the last subsection. Nevertheless, wc wish to emphasize again that the strong interaction between heavy quarks 
seen in this expansion is indeed also a feature of the underlying functional equation. 

VI. THEORETICAL ASPECTS OF CONFINEMENT 

A. Absence of long range forces 

It is sometimes argued that a strong, confining, quark-quark scattering kernel based on non-perturbative one-gluon 
exchange is not possible because it gives rise to Van der Waals forces between color singlet hadrons. Confinement, 
however, implies the absence of these forces. In nature this corresponds to nuclear forces described by an exponentially 
suppressed potential of the form 

V n (R) oc e - m - R , (68) 

at medium to large distances. Prominent long range forces also are not observed in quenched lattice simulations 
[75| . We therefore have to address the question whether the strong quark-quark scattering kernel found in this work 
necessarily implies power-law suppressed R~ aa Van-der- Waals forces between hadrons. As will be argued in the 
following, this is not the case. 

Let us begin with simple examples, building up to the more general case. First note that a single gluon exchange 
between color singlet hadrons is forbidden by color algebra. Next consider a two-gluon exchange between the two 
singlets. A corresponding diagram is e.g. given by 
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(69) 

where the two gluons combine to a color singlet exchange. Such diagrams arc nonvanishing and indeed expected to 
drive pomeron exchange [76|. However, they could only lead to Van der Waals forces in the absence of a mass gap, 
implying a massless glueball in contradiction with lattice gauge theory. Consequently such diagrams should give rise 
to a Yukawa potential Vn ~ e - M oR m £ ne s t a ^j c limit , which is suppressed at long distances R by the glueball mass 
M G . 

One can naturally employ a strong one-gluon exchange mechanism if simultaneously quarks are exchanged to 
balance color, namely 
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where we have explicitly shown the Bethe-Salpeter amplitudes coupling quarks to the physical hadron singlets. If we 
evaluate this diagram in the center of mass frame the two quark and gluon vertices and all propagators are evaluated 
in the small momentum scaling region. As a consequence, the infrared counting would give a power of (p 2 ) -2 from 
the exchange kernel, which is countered by powers of (p 2 ) 4 + 4BS from the two momentum integrations and the Bethe- 
Salpeter amplitudes. However, these amplitudes are finite in the infrared (BS = 0), since the Bethe-Salpeter equation 




displays a cancellation of the type (p 2 ) 2-2 between powers of momenta from the enhanced quark-quark scattering 
kernel and the internal momentum integration. Finite Bethe-Salpeter amplitudes are also required from physical 
considerations within the context of core nucleon-nucleon repulsion [77j . We therefore conclude that also the amplitude 
(|70| does not lead to colored Van-der-Waals forces. 

This above analysis is quite general. For the purpose of analyzing long-range forces we can limit ourselves to 
color singlet exchanges as the only non-vanishing amplitudes. However, in these diagrams the absence of long-range 
residual forces can be guaranteed by the sufficient condition that all color singlet, full scattering kernels in the theory 
are gapped, that is, they have no pole with zero mass. This condition is natural in that no massless particle is present 
in the experimental spectrum. Of course its validity remains to be shown in full Quantum Chromodynamics. Our 
point here is only that enhanced one-gluon exchange potentials do not necessarily generate Van der Waals forces 
between singlet hadrons. The detailed behavior of the color singlet kernels that appear in diagrams like ([55]) and ([70]) 
remain to be investigated in future work. 

Finally, it is worth examining the original Hamiltonian argument for long ranged colored forces given in [78| . One 
divides the Hilbcrt space into two sectors, P and Q, where P projects over states in which two interacting quark 
and gluon clusters are in a color singlet configuration and Q is the part in which the clusters are in non-singlet color 
configurations (the total state still being a singlet). Then the Hamiltonian can be exactly reduced to the color-singlet 
P sector as 

PH eS (E)P = PHP + PHQ E _ l QHQ QHP . (71) 

The presence of an infrared enhanced confining potential aR in the unconfined-color sector QHQ appears to yield 
a power-suppressed PH c gP oc 1/R Van der Waals potential in the singlet-singlet (confined) sector. However, the 
correct procedure to obtain non-relativistic inter-hadron interactions from the Landau gauge Green's functions is to 
perform a non-relativistic reduction not of the quark-quark but of the hadron-hadron scattering kernel. Since the 
kernel is gapped in the infrared, such a procedure would lead to Yukawa interactions 

1 _> 1 _> e -MR 

k 2 + M 2 k 2 +il/ 2 

It is certainly not correct to perform a non-relativistic reduction of a confining quark-quark scattering kernel and then 
attempt to construct a hadron-hadron interaction Hamiltonian from it. Thus the formal arguments of Gavela et al. 
do not apply to our approach. 



B. N-ality dependence of the string tension 

There are good arguments from both continuum and lattice QCD that the string tension a between static (anti-) 
quarks follows scaling laws that are dependent on the distance R between the quark and antiquark and the representa- 
tion of the gauge group. At intermediate distances R the string tension between quarks in the representation r follow 
the so called Casimir scaling, i.e. a r — C r /Cpo~F, where the subscript F denotes the fundamental representation 
and Ct,Cf the respective Casimir operators. At asymptotically large distances, however, the scaling law changes to 
°~r = f(k)ap, where f(k) is a function of the N-ality k of the representation r (see [ill 80j for reviews). 

Casimir scaling is exactly derivable from one-gluon-exchange type forces, thus the first term of the skeleton expansion 
in fig. [20] agrees with this behavior. All other diagrams involve quark propagators. From our discussions in sections 
IIIIHI and II V CI we recall that for really heavy quarks we find the perturbativc mass dependence of the quark-gluon 
vertex and the propagators implying that all these diagrams are suppressed by powers of the quark mass M. Thus 
approaching the static limit we recover Casimir scaling in our quark-antiquark potential, as we should [8lj . 
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Next we discuss the N-ality dependence of the asymptotic string tension. The N-ality of a given representation is 
the number of boxes in its Young tableau, modulo N. For SU(3) one has k = 1 for the fundamental representation 
(one box for the quark) and k = for the adjoint representation (3(mod 3)=0) for the gluon or octet fermion. In 
string- type models inspired in lattice gauge theory such as the Flux tube model [82l . l83l . Hj] , the string is postulated 
to snap if the hadron may decay to two color-singlet clusters. The energy stored in the flux tube has to exceed the 
threshold for the hadron-hadron combined energy and thus this happens only at large distances. It then follows from 
color algebra that the string between two fundamental sources in a triplet /anti-triplet representation can never snap 
in quenched approximation, since sea quarks are absent. Only gluons in the adjoint representation can be produced 
from the vacuum and these cannot combine with fundamental quarks to yield a color singlet hadron. However the 
octet-octet string can break, via two-gluon production, one of them combining with each source to make a color 
singlet. Thus, it is apparent that the string tension indeed scales with the N-ality of the representation. 

It is instructive to discuss how this works in the Green's functions approach to Landau gauge QCD. Certainly, 
the quark-quark scattering kernel with one gluon exchange cannot be responsible for this type of behavior, since it 
also confines fcrmions in the adjoint representation, leading to Casimir scaling. However this is only a four-point 
Green's function and by far not the only contribution to the computation of the string tension, say, via the Wilson 
loop. The following two graphs represent the first skeleton diagram of the 6-point Green's function that make N-ality 
scaling conceivable. The left diagram corresponds to fcrmions in the adjoint representation, the right diagram in the 
fundamental representation. 




If we send all scales simultaneously to zero as P — > 0, both diagrams are strongly infrared enhanced, with the help 
of our counting we establish a power law (P 2 )~ 3 ~ 2k . However in this case, all colored particles are simultaneously 
being separated, and we would still expect confinement to be at work. If now we separate the two clusters, keeping a 
finite momentum between the components of each cluster, most vertices and propagators have finite momenta running 
through, and only the middle gluon rung contributes to an infrared exponent of (p 2 ) _1+2K , therefore being suppressed. 

It is obvious that a full evaluation of the string tension from the Wilson loop requires the interplay of many such 
higher Green's functions, with complicated color structure. The interplay of these to yield the correct N-ality is a 
complete other project in itself. However it is obvious from these comments that demanding N-ality to arise from our 
quark-quark scattering kernel alone is not appropriate. It should (as it correctly does) yield Casimir scaling. 

VII. PHENOMENOLOGICAL ASPECTS OF CONFINEMENT 

A. Kernel non-relativistic reduction 

We would now like to examine the limit of our quark-quark scattering kernel when the quark mass is large and see 
how close we can get to a Schroedinger description. Before we proceed, let us remark that within full QCD a description 
of the quark-antiquark interaction in terms of a time independent potential is not entirely adequate (85l . l86l . l87j . 9 
The time-independent Schroedinger equation following from Coulomb-gauge QCD would be formulated in a Fock 
space with explicit transverse gluons. However, considering the strict non-relativistic reduction of the quark-quark 
scattering kernel in Landau gauge, i.e. seeking for a Schroedinger equation for a two-body problem, we need to insist 
in the elimination of all gluonic degrees of freedom. This then leads to a strongly energy-dependent kernel, which is 
equivalent to a time-dependent potential in coordinate space. 



One argument is: Within such a description the extension of a heavy meson at fixed radial quantum number n r would become vanishingly 
small in the static limit, i.e. < r > nr — > 0, and see exclusively the Coulomb part of the potential. This part being scaleless all Balmer 
splittings m nr+1 — m nr would scale with the quark mass M . However, the glueball masses wlq would not. So the states neighboring 
any one n r heavy quarkonium state would not be the corresponding n r ±\ state, but instead a state with mass m„ r + vtlq. We thank 
Gunnar Bali for this observation. 
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Nevertheless it is instructive to go through the exercise, since there may be a window of quark masses where a 
potential could be of use (while certainly important effects from virtual states with transverse gluons are expected 
too). This allows us to make contact with time- honored phcnomcnological and lattice computations in the static 
limit. 



B. Salpeter equation 



The main problem in deriving a Salpeter equation for a bound meson state are energy singularities in the loop 
momentum fc of the Bethe-Salpeter equation. In QED, ladder and cross-ladder diagrams in the Bethe-Salpeter kernel 
conspire to cancel these singularities, as detailed in appendix [j] However, due to factors from the color algebra such 
cancellations are unlikely to also occur in QCD. If one insists in maintaining the simplicity of the ladder approximation 
one therefore faces singularities in fco from both the quark pr opagato rs and the kernel (gluon propagators and vertices) . 
In non-relativistic reductions it has been common practice pa. l89l. fooL l9lj , to ignore the fco dependence of the kernel 
and Bethe-Salpeter wavefunctions and perform the pole integrals by putting only the intermediate quarks on-shell. 
The Bcthc-Salpctcr equation is thus inconsistently evaluated with ko = for the kernel and fco ^ for the propagating 
quarks. This procedure is justified empirically by demanding finite support for the BS-amplitude when the quarks 
are off-shell. 

The resulting Salpeter equation has many nice features such as color confinement (by sending colored objects to 
infinite mass outside the spectrum) in the presence of infrared enhancement kernel. It also respects chiral symm etry 
and it can be obtained as a collective excitation (Random Phase Approximation) in a Hamiltonian framework [92l. l93j| . 

In the following we go through the derivation of the Salpeter equation and discuss how it can be used to gen erate 
a static potential from our relativistic quark-quark interaction kernel. Here we follow Bicudo and Ribeiro [8!|; see 
also Itzykson and Zuber, [7l| . for an analogous derivation without propagator dressing functions), to define directly 
a positive (negative) Salpeter function from the Bethe-Salpeter amplitudes to emit (absorb) a quark-antiquark pair 
from the meson, evaluated at zero energy and integrated over the propagators 
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(72) 



(and analogously for ). The Bethe-Salpeter equations then reduce to the coupled Salpeter equations for the two 
amplitudes (j)^, the first of which becomes, now including the spin index, 
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(73) 



where the interaction part can be represented graphically with an instantaneous Bethe-Goldstone diagram (note 
fermions and anti-fermions are distinguished, anti-fermion always the bottom line) 




The integral over the propagators in the first line yields 

1 

M — Eg — Eg 
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and puts the quarks "on-shcll". One should note however that in the presence of a strong kernel, E q = 
\J A 2 (q)q 2 + B 2 (q) is actually divergent [H, [94[ (when the same approximations are applied to the propagator 
Dyson-Schwinger equation) and these divergences in the self-energies cancel only for color singlet states. 

To take the strict static limit one neglects « (j) + (this approximation is equivalent to Tamm-Dancoff and 
breaks chiral symmetry). Then Fourier transforming to coordinate space by means of 

v{r) = 1 w? v{q)elv (74) 

and 

Ei(r) = J (03^ r (^(g) - ^(0)) (75) 

one finally obtains a Schrocdinger-type equation (the infrared cancellation for color singlets is written explicitly) 

Aty(r) = J d 3 r' [E q (r - r') + E 9 (t - r')] + WW - E q (q = 0) - E q (q = 0)ty(r) . (76) 

The potential V(r) in this equation will be considered further in the next subsection. 

The insight one gains from the Salpeter construction so devised is that the static potential obtained from a Born- 
approximation is correct (modulo the infrared subtraction as needed for bound state equations with self-energies 
included) 10 . Thus one can test the quark-quark kernel constructed in Landau gauge QCD by evaluating the static 
potential in the Salpeter approximation. However, it is difficult to estimate the systematic error of this procedure. 
One way to estimate this is to cross-check with the Gross approximation, detailed in appendix IJ 21 Relativistic 
corrections to the potential of 0(v 2 ) could be investigated by appropriately expanding the spinors contracting the 
kernel and employing a Fierz rearrangement. These could also be evaluated in the Salpeter approach and compared 
to results from the lattice In general, the systematic error of the Salpeter equation is controllable in 1/M if one 
can independently establish that the BS wave amplitude has limited support in the quark off-shellness 11 . 



C. Salpeter string tension 

We are now using the results of the last subsection to relate phenomenology with the fundamental building blocks 
of our quark-quark scattering kernel, i.e. the quark-gluon vertex and the gluon propagator. To this end we will take 
the limit of our kernel with the quark-legs going on-shell, yielding qo = for the exchanged energy (ignoring possible 
singularities of the amputated kernel). The resulting non-relativistic potential is the a ppr opriate one to be used in 

3(1 1 C /~^T/2 



the context of the Salpeter (RPA) equation and is known to be about a = |0.18 GeV [92j, |95j. 

We denote the momenta of the incoming quarks by p\, k\ and the outgoing p2, The momentum transferred 
by the gluon rung is q. We choose the center of mass frame, where pi = — ki. We let W be the total energy and 
Pi = W/2 + £7/2, ki = W/2 — E/2 the energy of the two incoming quarks (E' for the outgoing). From the on-shell 
conditions for the quarks one finds that E = E' = in the center of mass frame. The gluon momentum is given by 
g=(0,q). 

Pl = (f'Pl) P2 = (f,P 2 ) 



g=(0,q) 

fci = (f,-pi) fc 2 =(f,-p 2 ) 



(77) 



Moreover one learns that in deriving cq. H76H a collective mode, namely the back-propagating Salpeter wavefunction, has been neglected. 
Taking this mode into account gives rise to additional energy-dependent hypcrfine-terms proportional to 1/M 2 in the non-relativistic 
Schrocdingcr equation. 

We thank A. Szczcpaniak for calling our attention to this point. 
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The one-rung skeleton amputated kernel (|77|) then becomes 



12 

K IR = g 2 C F rf(pi,p 2 )A J (M 2 ,q 2 ) D^(q) T^(k u fc 2 )A i (M 2 , q 2 ) (78) 



'•7 = 



where the Casimir CV = "^jy 1 stems from the color trace and the quark-gluon vertex is denoted by = X)l=i Aj. 

Next we contract each vertex of the kernel with external spinor wavefunctions (u^T^u) / Z 2 (should the kernel be 
part of a loop, the propagator attached to the right vertices needs to be decomposed as S = uu^/(D + ) + vv^ /(D_), 
and we keep the positive energy pole) . Note that the spinors need to be multiplied with the inverse of the wavefunction 
renormalization constant Z 2 in order to make the construction renormalization-group invariant. With the intent of 
obtaining the central potential we average over initial and sum over final spins. This means that at each vertex we 

have i XaiAo (^iiXzU^T^ jZ 2 . In the extreme non-rclativistic limit, the 1/M lower components vanish. In terms 
of Pauli bispinors we are then left with 



2 



0000. 



Herein, we have made a choice of quantization axis for the spin so that rotational invariancc is no more manifest. 
If necessary this can be avoided by employing a generic vector n. Plugging this spinor combination into eq. ([75]) 
and tracing each of the two vertices (with twelve tensors each) is conveniently carried out with the help of a symbol 
manipulation program such as FORM [96| . We furthermore define pi = p and perform the Wick-rotation to Euclidean 
space. 

Employing momentum conservation, k 2 = k± + q, p 2 = pi — q, only terms with the four "longitudinal" Xi pieces 
remain, with all transverse tensors Tj contributing to subleading parts of the potential. For the A^ we may write 

X i (M 2 ,M 2 ,q 2 )=A i (q 2 )- 1/2 - K , 
using our results from previous sections. After some algebra we are left with 

Kin = Ul + 2A 1 W(A 2 - A 3 ) + 4(WA 2 - A3) 2 Up • q) 2 - p 2 - AA\ ((q x p) 2 - (p 2 - pi))) (80) 



where the factor Z is the coefficient of the power law Z(p 2 ) = Z x (p 2 ) 2K of the gluon propagator in the infrared. 
Note that the total energy W ~ M and p cx 1/M such that the last term is suppressed for large quark masses and 
the frame dependence disappears from the central potential. To extract the central, static potential in the center of 
mass frame, we keep the first three terms and take W = 2M. Comparing ([50)) to the Fourier transform of a linear 
potential ar — > Sira/q , we read off 



a s Z 



(Aj + 4MA 1 (2MA 2 - A 3 ) - 4M 2 (2MA 2 - A3) 2 ) (81) 
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and a s = g 2 /Air. 

This string tension has to match the one extracted from the quarkonium static potential, a exp = (0A2AGeV) 
Z 2 is very close to 1 and Z ~ 2Q(GeV)~ AK can be read- off the solution for the vacuum polarization DSE, with a 
renormalization point fi ~ ISGeV, such that a s (fi) ~ 0.5. Ignoring A 2 and A3, we then have a phenomcnological 
estimate A\ ~ 0.16 GeV 1+2K where the energy units are GeV but the precise mass dimension depends on k, and is 
here fixed from the gluon propagator Z to k = 0.595. 

Of course, this number is just an order-of-magnitude guess, since the string tension obtained from the quark-quark 
2PI needs to be obtained by summing a skeleton series as previously discussed in length. The second term in the 
skeleton is the crossed-rung, one-loop diagram that causes a correction to the one-gluon string tension 



\ 6 M 2 J 



(82) 



where c is a pure number containing integration factorials, phase space and angular factor. We explicitly show the color 
factor, and M is a quark mass scale showing convergence of the skeleton expansion in its presence. The corrections 
to the string tension can be organized as 

ac^(l- Cl ^ + ...) (83) 
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To conclude this subsection, we note in passing that in the light of our calculation the longstanding discussion on 
whether confinement is of a vector 77 or of a scalar nature seems oversimplified. Each of the quark-gluon vertices in 
our expansion of the quark-quark scattering kernel contains a full set of twelve Dirac tensors with vector, scalar and 
tensor structures. After contraction with the gluon propagator one obtains a large number of terms and non-trivial 
cancellations can arise. The detailed structure of the resulting quark-anti-quark interaction needs to be investigated 
further and is certainly of great interest due to recent improved lattice data [97[ . 



D. Mass dependence of the string tension 

In this subsection we consider the possible mass dependence on the string tension from a phcnomcnological point 
of view. To this end it seems sufficient to consider the non-relativistic qq Schroedinger equation, that contains some 
of the features exploited in more elaborate schemes as for example pNRQCD [98[ . The reduced radial wave function 
u(r) satisfies 



h 2 

—u"(r) 
2p v ' 
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l(l + l)K 
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u(r) = (84) 



with a power-law potential V(r) = \r u and the reduced quark mass p = Mq/2. In particular, v = 1 and A = — a 
corresponds to the linear potential, v = — 1 and A = a s to the (attractive) Coulomb potential, and the limit v — ► can 
be used to represent a logarithmic confining potential. The equation can be transformed to a scaleless form [9{J Il00j | 
by performing the change of variables (r, E) — > (p, e): 

Therefore all the mass and parameter scaling of the spectrum and momenta (inverse length) needs to reside in 
cq. (|85[) . For the Coulomb potential v = — 1 one obtains 

p oc pa s , E oc fia s , (86) 

stating that all splittings grow with the quark mass at large Mq since the (low-lying) states become smaller and see 
mostly the short-range Coulomb potential of pQCD. A linear potential, however, causes splittings slowly decreasing 
with the quark mass and growing with the string tension 

p oc (pa) 1 ' 3 , E oc pT^a^ . (87) 

In fig. (|23p we show splittings between the 2s and 3s quarkonia for T, ip, and pir systems. Ignoring s — d wave 
mixing a typical assignment for these states are T(10355), T(10023), 004040), -0(3686), 0(2170), 0(1680). There is 
no candidate for I = cj(3s) so we give instead the I — 1 n, p splitting using p(1900), p(1465), 7r(1800), 7r(1300). 
However, one has to bear in mind that large errors in the mass determination and the general failure of a non- 
relativistic quark model for light quarks make these splittings less useful. The thick line in our plot represents a 
scaling based on eq. (|85() where a mass independent string tension is assumed, i.e. v = 1. 

With respect to these data we now discuss whether a strong mass-dependence of the string tension is possible. The 
and "0 states fall nicely on the p -1 ' 3 dependenc e sug gested by a linear-potential with a mass-independent string 
tension (thick line). However, as also noted e.g. in [105| . charmonium and bottomonium have a congruent spectrum, 
i.e. the splittings arc almost the same in both systems. This suggests either (i) a logarithmic potential, (ii) a soft 
mass-dcpcndcnce of the string tension of a linear potential of the type a oc U 'Mq, (iii) the beginning of a transition 
to the Coulomb regime around the T. 

On the other hand, our analytical analysis of the mass dependence of the quark-gluon vertex in subsection IIII HI 
yielded three possibilities: either the vertex-dressing is proportional to Mq which results in a oc Mq, or it is pro- 
portional to 1/Mq, yielding a oc Mq 2 , or it is independent of the quark mass which gives a oc Mq. The first case 
together with eq. (|87|) would give a Coulomb-like behavior E oc Mq that is clearly ruled-out by the data. The 
decreasing self-consistent solution a oc 1/Mq would give E oc Mq 5 ^ 3 that is too strong a decrease for the data. 

On the other hand, the third solution has a quark mass independent string tension, scaling as the square of the 
QCD scale, just like in pure Yang-Mills theories [l06l |. This behavior could explain the experimental splittings if the 
Coulomb potential has a large effect in the bottomonium spectrum. It would certainly be desirable to confirm this 
behavior in a lattice computation of the dependence of the string tension with the quark mass. 
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FIG. 23: Left diagram: The dependence of the splitting E(3s) — E(2s) on the quark mass. The experimental assignments for 
uu-, ss-, cc-, and bfe-states are educated guesses based on the non-relativistic quark model. The thick line represents a scaling 
based on eq, (185|l where a mass independent string tension is assumed, i.e. v = 1. 

Right table: Charmonium and bott omonium spe c tra (quoted are the splittings to the J/i/)(3097) and T(9460)). Data from 
[Hal, with new states reported from [lOll . Hoi [io3 . Ho3 |. 



Moreover, many of the experimental splittings given in the table in fig. [23] agree with a quark mass-independent 
scaling, including many of the exciting new states. Of course, in the absence of coupled open decay channels one would 
like to make the comparison for the highest possible states, since the linear potential would then clearly dominate, 
but the open channels make the loss of precision high in the spectrum comparable to the loss of precision due to the 
Coulomb potential for low-lying states. Also, some of the states could be hidden exotics in nature, that is, composed 
of wavefunctions qq q, qqqq an d one needs careful on-resonance data to be able to apply the Franck-Condon principle 
to distinguish them [l07l . [l08j . 



VIII. SUMMARY AND CONCLUSIONS 



In this work we performed a thorough investigation of dynamical chiral symmetry breaking in the quark-gluon inter- 
action of quenched Landau gauge QCD. We considered the coupled dynamical system of Dyson-Schwinger equations 
for the quark-propagator and the quark-gluon vertex and determined selfconsistent solutions both analytically and 
numerically. We identified the infrared leading diagrams in this equation and found these to agree with the leading 
ones in an l/iV c counting scheme. This result is alternative to a number of treatments of the vertex-DSE (see e.g. 
[TBI [l6l, l4(il. [4l| . that picked out the subleading, Abclian parts of the vertex-DSE. The most important results from 
our selfconsistent treatment are 

• The qqg vertex uniform singularities have, in a milder form, been proposed in the past, in Ball-Chiu-type 
constructions based on the Slavnov- Taylor Identities. We now find an overall (all external scales go to zero) as 
well as a collincar (the external gluon momentum goes to zero) infrared divergence of the vertex scaling with 
(p 2 ) -1 / 2-K , provided we take into account dynamical chiral symmetry breaking. This includes corresponding 
chirally non-symmetric tensor structures in the quark propagator and the quark-gluon vertex, as e.g. Dirac 
scalars. In our numerical analysis we find a scalar part of the vertex which is of the same strength in the 
infrared as the conventional 7^-part. 

• This infrared divergence implies an infrared divergent running coupling (defined from the quark-gluon vertex) in 
contrast to the infrared finite fixed point behaviors of the couplings in the Yang-Mills sector of the theory. Our 
findings agree with the old notion of infrared slavery [l[, although the mechanism at work is different. In early 
works it was believed that infrared slavery is caused by an infrared divergent gluon propagator. In contrast we 
find an infrared vanishing (or finite) gluon, that is accompanied by an infrared divergent quark-gluon vertex. 
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Both quantities together drive the coupling to infinity in the zero momentum limit 12 . 

• The collinear infrared divergence also triggers an infrared divergent one-gluon exchange kernel between two 
quarks, leading to a linear rising static quark potential which implies quark confinement. 

• If chiral symmetry is restored by hand, i.e. if the chiral symmetry breaking tensor structures in the quark 
propagator and the quark-gluon vertex are arbitrarily set to zero, we find a different selfconsistent solution of 
the quark propagator and vertex-DSEs: the vertex then has an infrared singularity of order {p 2 )~ K which is 
reduced in strength compared to the chirally broken case. 

• This reduced infrared strength in the chirally symmetric case leads an infrared constant coupling from the 
quark-gluon vertex (i.e. similar to the couplings in the Yang-Mills sector) and to a Coulomb-type static quark- 
antiquark potential, which is not confining. We have therefore uncovered a mechanism which directly links 
confinement to dynamical chiral symmetry breaking and vice versa. 

These are the rigorous results reported in this paper. We have also, in a less rigorous but nevertheless meaningful way, 
discussed some aspects and consequences of our findings. In particular we gave a prescription to construct a confining 
and chiral symmetry preserving Bcthe-Salpeter kernel, which includes the infrared divergencies found in our analysis. 
We discussed the mass dependence of the resulting string tension in our approach and found a surprising behavior 
for the quark mass functions: taking all tensor structures of the quark-gluon vertex into account it seems that the 
infrared results for the quark mass M(p 2 = 0) are much less sensitive to the current quark mass than anticipated 
previously. Whereas there is a change of O(10 3 ) in the current quark mass from the up- to the bottom quark, we 
only find changes of a factor 2-3 in M(p 2 = 0) (somewhat larger changes could be expected in different truncations). 
As a result we find a string tension which has no simple scaling law wrt. to the quark mass. This result seems to be 
consistent with phenomenological expectations from the splittings in quarkonia. 

Finally we discussed some aspects of the confining nature of our result. We presented arguments in favor of the 
absence of long-ranged van-der-Waals forces even in the presence of a strong singularity, such as the one suggested, 
in the quark-gluon interaction. In a nutshell, the reason here is the absence of color-non-singlet exchanges between 
color-singlet objects. The remaining color-singlet exchanges then always are gapped such that the interaction vanishes 
at low momenta. This issue deserves further studies in the future. In addition we shortly discussed the Casimir scaling 
property of one-gluon exchange and potential candidates of Green's functions which should lead to the N-ality behavior 
of the string tension at large distances in unquenched QCD. 

Finally, it is hard to se e ho w the picture of confinement via a static linear potential may be also adequate for light 
quarks (see for example jl09j ]). In this respect it is reassuring that we find positivity violations for the light quark 
propagator, which in turn is a sufficient condition for the confinement of light quarks. At the current technical level 
of treating the DSE of the quark-gluon vertex this statement is truncation dependent. It remains to be seen, whether 
this finding can be corroborated in a truncation independent way. 

In this work we gave a qualitative picture of what could be a confinement mechanism for quarks in Landau gauge 
QCD. Although we tried to make our arguments as rigorous as possible, much work remains to be done. In particular 
it remains to be shown that the quantitative aspects of the picture match expectations from phenomenology. Here it 
seems desirable to develop a calculational scheme to determine the quark-antiquark interaction in quantitative detail 
and confront the results with the experimental findings as e.g. the spectra of charmonia and bottomonia. 
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APPENDIX A: STRUCTURE OF THE QUARK-GLUON VERTEX 



The quark-gluon vertex is in general a function of 2 independent four momenta as well as its Lorentz and color 
indices. The tensor decomposition expresses it in a finite basis of the internal spaces given in Sec. IA 11 This also 
reduces the momentum dependence of the form factors to 3 independent momentum squares. We will study several 
different kinematics defined in Sec. IA 21 



1. Tensor basis 



In Landau gauge the quark gluon vertex has twelve different tensor structures conveniently written as 

4 8 



T^ = J2^ + J2nTt , (Al) 

i=l i=l 

where the Ti denote parts transverse to the gluon momentum. The general basis for the quark-gluon vertex we employ 
here has been given by Ball and Chiu [24[ . It reads 

Lt = r (A2) 

l» = -i(jp l+P2 y 

T l = i(PlP2-P3-P2Pl-Ps) 

T 2 = iPlP2-P3 - P2Pl-P3)(jPl+ fa)) 

n = fap%-ph» 

U = ^(P3)u 

n = (A-pt)r+(pi+P2r h 

7? = \{ P \-Pl)[{fa+ fa)l»-{pi+P2Y]-i{pi+P2YaxvP^p l i 

Tg = -^a Xv P2P\- faPi + faP2 

with P 3 = P2 — Pi and = 5(7^7,, — 7^7^)- This basis is similar to that of [ll|, differing only in the momentum 
labels. The structures L\, L 2 , T 2 , T 3 and T 6 each involve one Dirac matrix and therefore form the Dirac-vector part 
of the vertex. The structures L3 and T\ form the scalar part of the vertex, whereas the remaining structures L4, 
X4, T5 and Tg have tensorial character. For the purpose of our power counting analysis as given in section Dl the 
tensors L4, T4, T5 with an even number of Dirac matrices act similar to the scalar ones, whereas the tensor Tg with 
an odd number of Dirac matrices acts similar to the vector part. 



2. On the kinematics of the quark-gluon vertex 



Since the quark-gluon vertex Dirac amplitudes Xi and n, defined in eq. ((2]), are functions of three Euclidean scalars 
Pqi P% P 2 g i it is customary to plot them in slices of this three-dimensional space, that we denote "kinematic sections". 
Given the various choices in the literature, we find it convenient to collect here definitions for the most used sections. 

First, there are important kinematic sections that are "degenerate", in the sense that there is only one independent 
vector p M in four-dimensional Euclidean space instead of the usual two vectors p^, p 1 ^ (the gluon momentum being in 
the plane spanned by them). 

A popular degenerate section is the "soft gluon" section, characterized by a gluon of exactly zero momentum 
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P /\p 

A similar section is the "soft quark" section where one quark momentum vanishes. This corresponds in the chiral 
limit to an on-shcll quark in Minkowski space: 

P 

o/\ P 

Next let us look at two interesting non-degenerate sections. We introduce the so called "1-2-3" or "totally 
asymmetric" section by 




imposing the restrictions (in terms of p\ = p 2 ), p\ = 2p 2 , p 2 = 3p 2 that fix the relative angle of pi, p2 and the quotient 
of their Euclidean moduli. 

The other section is the so called "symmetric" section given by p\ = p\ = p\ = p 2 , in which all scales flowing into 
the vertex have equal value (but p\ and pi are not parallel) . 

These four sections are in principle useful to study the infrared behavior of the vertex, that is, the infrared exponents. 
To expose the angular dependence of the vertex dressing functions, it is convenient to define a "fixed scale" section 
in which we vary the angle between p\ and p2 but fix p\ = p\ = p 2 ^ 0. Thus we can plot \i(x) in terms of the polar 
cosine 

^ = Pi -P2 p\ 

x M H 2 P 2 ■ 

APPENDIX B: SECOND DSE FORMULATION FOR THE QUARK-GLUON VERTEX 

In section irm in the main part of this work we gave a solution for the infrared behavior of the quark-gluon vertex in 
terms of power laws. There it turned out that the vertex has both on overall infrared singularity of degree (p 2 ) _1 / 2_K 
when all external momenta are proportional to one scale p which goes to zero and also a collinear singularity (p 2 )^ 1 ^ 2 ^ 1 ^ 
when only the gluon momentum p g is small compared to Aq C d- 

This solution has been derived from a version of the vertex-DSE, where the bare vertex is always the one attached 
to the external gluon leg. This gives the equation shown in fig. [3] in the main body of this work, which includes all 
possible bare vertices from the Lagrangian of the theory. We will call this version of the DSE 's-channel version' in 
the following. This solution also straightforwardly solves the corresponding equation of motion for the vertex from 
an 3PI effective action, given in fig. [5j However, things are slightly more complicated for the second or 't-channel' 
version of the vertex-DSE. This version has a bare vertex attached to an external quark line and is given in fig. [3J 

The problem here is that the simple triangle diagrams involving only two- and three-point functions do not reproduce 
the infrared singularity found from the other two equations. Let us carry through this exercise for example for the 
diagram involving the fully dressed three-gluon vertex. Simple power counting then results in an overall divergence of 

(p2)-3„+l/a x (p 2 )2K -l x (p2)2«-l x (p2)-l/2-« x (p2) 2 _ {p 2 )0 ^ (m) 
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FIG. 24: Leading contributions in the skeleton expansion of the t-channel- version of the DSE for the quark-gluon vertex. 




FIG. 25: Kinematic configurations of the leading loop graphs that generate the strong soft-gluon divergence (a) and the uniform 
divergence (b) of the quark-gluon vertex. The labels s and h denote whether the momentum running through a line is soft or 
hard, respectively. In contrast to the s-channel-version of the DSE, here the soft-gluon singularity in graph (b) is necessary to 
induce the uniform divergence. 

where the first term stems from the three-gluon vertex, the next two terms from the gluon propagators followed by the 
scalar part of the remaining quark-gluon vertex and the loop integration, whereas the quark propagator is constant. 
To yield the previous IR divergence of the vector part of the vertex on the left hand side there is clearly infrared 
strength missing. 

The answer to this puzzle lies in the remaining diagram, that includes the 1PI gluon-quark four-point function. 
Since the ghost dynamics is the dominant dynamical contribution in the DSEs and there is no fundamental quark- 
ghost coupling, similarly to the case of the s-channel DSE discussed in detail in the main text, it is again necessary 
to perform a skeleton expansion to two-loop order as shown in fig. 1241 

The power counting analysis for this equation is complicated since not only soft modes of the order of the external 
momenta can dominate the different loop integrals but also hard modes of the order of the dynamically generated 
quark mass M appearing in the quark propagators 13 . In combination with the different possibilities to route any hard 
external momenta through the diagrams, the two-loop integrals can be decomposed into many different potentially IR 
sensitive regions (23j . The most singular contribution on the right hand side of the DSE dominates and will determine 
the scaling of the vertex on the left hand side. Naturally, Greens functions involving only hard momenta which appear 
in a given graph do not scale with the soft external momentum. Here, we will not give a complete power counting 
analysis analogous to the one for the first DSE presented in appendix [0 but rather identify the contributions that 
provide the strong singularities found in the s-channel DSE discussed in the main text. 

Let us start with the quark-gluon vertex in the soft-gluon limit. One particular contribution arising from a specific 



We note that a very similar integral appears in the DSE for the quark propagator, which is very well explored. These integrals are 
known to result in finite answers in the infrared for a wide range of possible dressings in the loops (see e.g. 
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kinematic region of the first two- loop graph in fig. [24] is given in fig. [25] (a). Here, the loop integral of the initial 
graph (left part) arising in the unexpanded DSE fig. 3] extends over hard momenta whereas only the loop integral 
(right part) arising from the skeleton expansion involves soft scales of the order of the external gluon momentum 14 . In 
case that this contribution dominates it yields for the anomalous exponent of the quark-gluon vertex in the soft-gluon 
limit the self-consistent result (counting only IR exponents) 

< = 2+ Q + <^)+ 2^ + 2(^-1) = 25*+ \+k => 8f g = - l --n (B2) 

Accordingly the ghost-enhanced three-gluon vertex induces a strong soft-gluon singularity in full analogy to the s- 
channel version discussed in the main text. However, this mechanism does not generate a uniform divergence in the 
absence of soft divergences. To see how the latter is generated self-consistcntly consider graph (b) of fig. [55] which 
presents one of the possible contributions. If it dominates it yields for the corresponding anomalous exponent 

5l = 2 + 25 u qg + 8f g + 8l g + {5 gl -l) + 28 q = 25 u qg + \ + K 8 u qg = -\-n (B3) 

We therefore conclude that the structure of the loop integrals in the last term of fig. [4] together with the collincar 
singularity in the external gluon leg provides exactly such strength in the infrared to reproduce our results for the 
quark-gluon vertex in the s-channel version of the vertex-DSE also for the t-channel version. 

It is interesting to note a few properties of this result: First, in contrast to the s-channcl DSE, here regions of 
the loop integral with hard momenta of the order of the quark mass dominate even in the uniform limit. Further, 
the uniform divergence in the t-channcl DSE is not directly induced by the three-gluon vertex but only indirectly 
once the strong soft-gluon divergence is present. Therefore, in the present case these two divergences are not totally 
independent but imply each other. Furthermore, in the t-channel version it is not necessary to distinguish between the 
different Dirac components of the vertex to recover the proper uniform scaling. The reason is that the corresponding 
graphs involve many fermionic Greens functions that involve hard momenta and do not scale with the soft external 
momentum, independent of their Dirac-structure. Correspondingly, in contrast to the case of the non-Abelian graph 
studied in the main text, there are no corresponding restrictions, here. 

Finally, from a practical point of view the s-channel version of the DSE lends itself much better to a numerical 
analysis, since, as explained in the main text, the IR dominant two-loop graph can be approximately replaced by 
a 1-loop graph with a dressed three-gluon vertex. In the quenched approximation this three-gluon vertex can be 
computed beforehand or, as done in the present analysis, simply be implemented by some ansatz. In the t-channcl 
version of the DSE this is not possible since the initial loop dresses the desired quark-gluon vertex itself and therefore 
the expensive two-loop integrals indeed would have to be analyzed numerically. 



Note that we are only interested in the scaling with respect to the soft external momentum p <C Aqcd and do not distinguish between 
any hard external momenta q < Aqcd and hard loop momenta q ~ M. 
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APPENDIX C: SEMI-PERTURB ATIVE ANALYSIS OF THE QUARK-GLUON VERTEX 

1. Non-Abelian diagram 

We start by evaluating the one-loop diagram involving the non-Abelian three gluon vertex 



*£ = (P2-Pl)" 



igT«r» A (pl^pl(p 2 - Pl ) 2 ) 




(CI) 



The color structure for this diagram is 



iN, 



2 ij 



where T°- is the generator in the fundamental representation, and the factor N c results from the contraction of the 
color matrices. The full expression reads 



-iN c 



(So-eifa ~ kOju + <W(&2 - h)a + S atJ .(k3 - ki) e ) I 8 ua lv . n la ) 



(C2) 



k\ 



with the renormalized coupling g 2 and the renormalization constants Z\p and Z% for the bare quark-gluon vertices 
and the three-gluon vertex. Dimensional analysis of eq. (|C2[) immediately shows a logarithmic divergence which 
appears only in Ai. All other pieces of the vertex, containing powers of either p^, are free of ultraviolet divergencies. 
Our renormalization prescription for Ai is given in subsection IIV Al 



2. Abelian vertex diagram 

Next we turn our attention to the "QED-likc" diagram 



igT^ A (pl 2 ,pi(p 2 - Pl ) 




ki = pi-q 



a%k% = (p 2 ~ Pl y 



-k 2 = P2-q 



(C3) 
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This has color structure 



2N r 13 



and reads 



Both the Abelian and the non-Abelian diagrams reproduce important properties of the vertex; as discussed in [Icl ] , 
charge conjugation invariance requires 

CYC' 1 = -7 M T CiT^C- 1 = iD" T 

where C is the charge conjugation operator and the superscript T denotes transposition. Referring to the WT basis 
(|A2[) above, we see that the A4, T4, and tq Dirac amplitudes should be odd under the interchange of p\ and P2, whereas 
the other nine should be even. This is indeed the case for both diagrams in our calculation. Also, the anti-unitarity of 
the time reversal operator means that all the Dirac amplitudes should be cither real or imaginary. With the explicit 
factors of i in cq. (P) and the conventions of (|A2|) all the A and r are strictly real. 



3. Input for gluon and quark propagators 

The dressing functions for the gluon Z and quark <r v , a s propagators are taken as input from reference |lldl | . We 

employ analytical fits to the numerical solution of the Schwingcr-Dyson equations from this work corresponding to a 

2 

choice a s = f^r = 0.2. The gluon propagator is well represented by: 

/ , x \ 1+25 

Z{x)=[^l) R 2 ( x ) (C5) 



and the ghost propagator by 



G{X)= {M) (C6) 
where the low momentum behavior is captured by the irrational function 

r{x) . ■m , :'tt) > , (CT) 

The parameters appearing in R(x), solution of the DSE are c = 1.269 and d = 2.105, and in the quenched approx- 
imation for the gluon and ghost equations 5(Nf = 0) = —9/44. The value for the exponent k ~ 0.595 in this fit is 
approximately known from the infrared analysis of the DSEs Q ■ 
Finally the running coupling is well reproduced by 

■M - —, M («) 



i»g( e '+« 1 ( ; &r+i.ite) 11 ) 



with ai = 1.106, a 2 = 2.324, bi = 0.004, b 2 = 3.169. Using a MOM scheme and fitting A to only the ultraviolet 
behavior, a value Aq C d = 0.714 GeV has been given in Ref. [9|. Note, that these fits are for faster computing and to 
offer a closed form: all results for the various two-point functions are solutions of their respective Schwinger-Dyson 
equations. These fits will be used for all numerical calculations presented in the course of this work. 

In the same way we can parameterize the quark propagator functions, a SiV (p 2 ) 1 and fit these to the numerical DSE 
solutions. Again, we want to make explicit the leading logarithmic corrections to Z^(p 2 ) and M(p 2 ), achieved with 

, 2. B chilal (p 2 ) C cqm a{p 2 + m 2 m 
a ^P > = 2 1 2 ' i L9 ) 

P + m S ing P + m si„g 

II ce(v 2 + m 2 :„„) 
av(P ) = P 2 + ™L S ( 1 " 2. + ° B ^ AP 
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with 

S C W(P ) - Q csb p2 + m 2 ^ + A2 + {p2 + m 2 ng + A2)2 ■ (CH) 

This parameterization ensures the correct asymptotic behavior for cr Sj „(p 2 ) and the quark functions M(p 2 ) and Z^(p 2 ) 
0, [53] in the chiral limit. It is important to note that we use this fit for the quark propagator in this appendix 
only! In the main text, where we back-feed the vertex into the quark-DSE, given in fig. [TJ the corresponding quark 
propagator is solved sclf-consistently and dynamically from the set of coupled DSEs for the quark propagator and the 
quark-gluon vertex. With this input for the propagators the loop corrections eq. (|Clj) and (|C3j) are straightforwardly 
evaluated and yield the result shown in fig. [7] in the main text. 



APPENDIX D: POWER COUNTING AND DIRAC STRUCTURE: AN EXPLICIT EXAMPLE 

In this appendix we demonstrate the fact that the power counting analysis depends only on the Dirac-scalar or vector 
nature of the individual tensors of the quark-gluon vertex explicitly in a special kinematic limit. Due to momentum 
conservation the vertex depends only on the two external quark momenta (pi)^ t , (p 2 )/j or three Lorentz invariants 
Pi?P 2 ?Pi ' P2- To analyze the infrared limit of this vertex in the presence of only one external scale p\ = p 2 <C A 2 CD , 
we can set (p 2 ) M = 2(pi) M ~: 2p^ without loss of generality (Note that the case (p 2 ) M = (pi)^ corresponds to the 
collinear limit where the gluon momentum goes to zero, while the other momenta stay finite. We treat this limit 
separately in subsection IIIIEI ) 

The choice (p 2 ) M = 2(pi)^ leaves us with only four possible tensor structures which can be denoted by 

4 

r^p)= t gJ2Up 2 )Gl (Dl) 

i=l 

with being Lorentz and Dirac scalar functions and 

Gj,=7„. Gl=p^ Gl=i>p^ G%=h^ (D2) 

where we have normalized the momentum, p^ — p^/ \fp 2 ~, to ease power counting. Note that G* and G^ have an 
odd number of 7-matriccs, whereas G 2 and G ^ have an even number. Therefore £2.4 7^ only if chiral symmetry is 
broken. 

To project out the respective dressing functions £i.. 4 (p 2 ) from the vertex-DSE we multiply the equation with 
Jfi,Pfii$Pfj,,fi"ffj, respectively and take the Dirac-trace. Linear combinations of these four basis elements then give the 
desired £i..4(p 2 ). 

For the internal quark-gluon vertices we can employ the expression (|D1[) with any internal momentum as argument: 
it contains all possible types of Dirac structures (vector, scalar, and tensor). Any more complicated dependence 
on external and internal momenta will generate the same powers of external momenta in £1. 4 after integration (for 
dimensional reasons). To determine the infrared exponents of the quark-gluon vertex we employ the ansatz 

Up 2 ) = xi ■ (p 2 )* , Up 2 ) = X2-(p 2 Y\ Up 2 ) = x 3 ■ (p 2 ) 03 , U(p 2 ) = xa ■ (p 2 )^ (D3) 

for the scaling of the dressing functions of the quark-gluon vertex with the momentum scale, where the pre-factors \i 
depend only on momentum ratios. 

Since we contracted the right hand side of the DSE with appropriate tensor structures, see above, we eventually have 
to deal with scalar integrals only, which depend on powers of internal and external momenta. The evaluation of such 
diagrams has been described in [561 ] - Most important for our purpose is the fact that all powers of internal momenta 
transform to powers of external momenta after integration. This is also clear from simple dimensional consideration, 
provided there is only one external scale present. In our case this is ensured by the condition p 2 <C A 2 CD < M 2 , which 
leaves only one scale p 2 in the deep infrared. The following analysis rests on the assumption that there is no second 
small scale present. In principle this assumption may not be true for heavy quarks, since then one has a potential 
additional small scale Aqcd/M 2 . We have investigated this possibility in our numerical calculations of the quark-gluon 
vertex and did NOT see any influence of such a quantity. The corresponding results of our quark mass study are 
presented in subsection IIV CI where we discuss this issue further. 
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A similar analysis can be done using all four dressing functions £i..4. The explicit integral equations are quite 
lengthy and add nothing to the power counting argument, so we omit them here. One then obtains: 

(pT-rnax \Z f — ,Z f — ,Z } — ,Z r 



M 2 ' 1 M 2 ' 7 M 2 11 M 2 

j-p2^1 + K+ft+/3 3 /p2\l+«+/3 2 +/3 4 /p2jl/2+K+/3i+j8a 



Z f~ 



M 2 ' M 2 ' M 

/•p2\l/2+K+0i+#4 ^2U/2+ K +ft+ft ^2^1/2+k+/3 3 +/34 -J 

Z/ M > Z * M ' Z/ M / (D4) 

(py^max | ^ ,Z f — 2 ,Z f — 2 ,Z f — 2 , 

(p2)l/2+K+2/3i ^2-)1/2+k+2/3 2 ^2)l/2+K+2/3 3 



' Z / 17 ' Z / 



J M ' 1 M ' J M 

£p2\l/2+K+2#4 ^2jl/2+K+0i+/3 3 fp2\l/2+K+^a+j84 "J 

Z ' M M ' Z/ M } (D5) 

, 2, P3 j 7 (p 2 ) 1+B+2ft 7 (p 2 ) 1+K+2fe 7 b 2 ) 1+K+2fe 7 (p 2 y+^ 

(p)f~max \Z f — 2 ,Z f — 2 ,Zf -= ,Z r 



M 2 1 1 M 2 ' 1 M 2 11 M 2 

^2^1+re+/3i+/3 3 (p2y+K+02+pi ^p2y/2+K+/3i+/3 2 



Z f 172 ' Z f 172 ' Z f 



M 2 ' M 2 ' M 

^2)1/2+k+/3i+/3 4 (p2)l/2+K+/3 2 +,33 ^2^1/2+k+/3 3 +/3 4 ~| 

Z/ M ' Z/ M ' Z/ M J (D6) 

(pY^max \Zj — 2 ,Z f — 2 ,Z f — 2 ,Z, — 2 , 

^2U/2+«+2/3i ( p 2\l/2+K+2/3 3 

Zf M ' Zf M ' 

( p 2\l/2+K+2/3 4 ^ p 2y/2+K+p 1 +p 3 fp2\l/2+K+l3 2 +04 ~\ 

Zf M ' Z/ M ' Z/ M (D7) 



which is solved by 

02 = -1/2- k, Pi, foe [-1/2 -k,-k] 

Pa € [-1/2 - k, 1/2 - k] . (D8) 

This leaves a range of possibilities with the common trait that the most infrared singular piece diverges as (p 2 ) _1 / 2 -« 
in the infrared. We have analyzed this analytic solution in numerical calculations, presented in section ITVl and found 
that only the case p2 = Pi = Ps = Pi = — 1/2 — k is realized. 



APPENDIX E: A SIMPLE MECHANICAL MODEL 



The simplest and well-known example of spontaneous symmetry breaking is given in fig. 1261 A parity-symmetric 
potential V{x) — x 2 — x 4 has two minima at x m = ±1/^2 and a classical bead at a minimum of the potential has to 
choose between left or right. Therefore the solution to the system in its ground state spontaneously breaks the parity 
symmetry, but the potential perceived by the bead remains parity symmetric. 



This example is analogous to the spontaneous chiral symmetry breaking in the Dyson-Schwinger equation for the 
propagator with a fixed quark-gluon vertex and gluon propagator, where a quark mass gap is generated as a function 
of the potential constants. 



In this paper we have considered a more general situation, and it is worth illustrating it with the following simple 
example. Consider now a bead of mass m deposited on top of an elastic band of elastic constant k and half-length 
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FIG. 26: A bead on a rigid ramp with two degenerate minima rolls down to either. The parity symmetry of the system is 
broken down spontaneously at the ground state. 




FIG. 27: A bead on an elastic band held in position at its middle point. The displacement of the bead (spontaneous parity 
symmetry violation by the "reduced" bead system in its ground state) is proportional to the stretch of the band (spontaneous 
parity symmetry violation by the potential that the bead sees), as measured for example by the area between the band and 
the horizontal. Both are coupled and it is not possible to understand one without the other. 



Iq. The half-band is held fixed at the two ends and the mid point. This situation is represented in fig. [57] The bead 
can of course stay at the middle of the band, with its weight compensated by the clamp holding the band by the 
middle, but this is a point of unstable equilibrium. A random perturbation will make the bead roll left or right. If 
the radius of the bead is small compared to lo, the band will adopt a triangle shape, with length li left of the bead 
and I2 right of the bead. The bead will be displaced from the horizontal a height y and from the center a distance x. 
The potential energy is now 



with the geometric constraints 



Minimizing the potential yields 



and, for small stretchings, 



k 2 
V(x, y) = -mgy + - ((h + h) - lo) 



x 2 +y 2 = l\; (l Q - x) 2 + y 2 = l\ 



dV(x,y) _lo 
dx ~ ~ 2 



9V(x,y) n 1 l2 ,,^1/3 



= y = ~ illmg/k) 



dy a 2 
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so that the bead in this minimum sees an effective potential that is approximately 

/ \ 1/3 

V(x) = -mg{-^-\ \h/2-x\ x > ; 

V{x) = x < 

reflecting the back-reaction of the mass of the bead on the stretchable band (if m = 0, V(x) = remains parity 
symmetric under x — > —x). 

The position of the bead is again analogous to the quark mass function in the propagator after spontaneous chiral 
symmetry breaking, and the slope of the potential (or equivalently the area between the band and the horizontal) 
analogous to a scalar piece in the vertex function, that also appears after spontaneous symmetry breaking. One 
cannot understand one without the other, and both are proportional. 

In a field theory of course, the number of degrees of freedom is infinite and not only two, and so are the coupled 
Schwinger-Dyson equations, all of which have to simultaneously display spontaneous symmetry breaking. 



APPENDIX F: POWER COUNTING WITH KINEMATIC SINGULARITIES 



In this appendix we present the full power counting analysis for the quark sector of quenched QCD, given by 
the coupled system of the quark propagator and the quark-gluon vertex DSE. In particular we take into account 
the possibility of kinematic singularities of the vertex and do not make assumptions on the IR limit of the quark 
propagator. The analysis of kinematic singularities follows the corresponding analysis for the gauge sector of the 
theory presented in [23 ]. 

Whereas the propagator is generally given by a single IR exponent S q , the quark gluon vertex has in addition to the 
uniform limit that all momenta vanish 6 qg also the possibility that only the gluon 8-* g or only the quark leg vanishes 
5 q g . As noted in subsection llll Dl in the main text it is crucial to distinguish between the scalar and vector components 
in the analysis. We distinguish them by additional lower case letters v for the vector and s for the scalar part, so 
that the full set of possibly distinct IR exponents is given by S qv , S qs , S qgv , 5 qgs , S^ gv , Sf gs , S qgv and S q gs . Where the 
exponents without a tensor specification are used, they stand for both Dirac structures so that terms with all possible 
combinations arise after decomposition. Finally, we define the canonical exponent of the quark propagator Xq which 
is given by \qs = and Xqv = 1/2 for the two cases, respectively. 

The distinction of the tensor components is already important for the quark propagator, where the different compo- 
nents feature a different scaling already at the tree level. In order to perform the power counting analysis for the DSE 
of the individual components of the quark propagator we first have to consider the scaling behavior of the appearing 
inverse of the dressed propagator which is given for the parameterization eq. ([3]) by 

S (p) = (ta v p 1 + cr s ) = i i>- -> i (p ) tf-(p) ■ (Fl) 

a(,p z + crj cr£p z + cr| v ' 

The DSE for the vector part of the quark propagator is given in fig. [28l where the labels V and S next to the 
propagator and vertices denote vector and scalar components, respectively. This DSE leads to an equation for the 



FIG. 28: Decomposed DSE for the vector part of the quark propagator. Vector and scalar components are denoted by V and 
S, respectively. 

corresponding IR exponent 

+ min ( 2 V + 1, 2 M = min Q, S^ gv + S qv + | + 2k, 5^ gs + S qs + 1 + 2k, S q qgv + ^,5« gs + ^ . (F2) 

Here it has been taken into account that due to a possible quark mass the loop integral can in addition to the IR 
regime also be dominated by hard modes reflected by the last two terms, where the 6% g arises due to the possibility 
of a soft quark divergence of the vertex and the 1/2 due to a further suppression arising from cancellations owing to 
gluon transversality. The corresponding DSE for the scalar part of the propagator is shown in fig. [55] which leads to 

S qs - min (25 qv + 1, 2S qs ) = min ( 0, 5 u qgs + 5 qv + | + 2k, 5^ gv + 5 qs + 1 + 2k, 8\ gs , 8 q qg A . (F3) 
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FIG. 29: Decomposed DSE for the scalar part of the quark propagator 



5 u qg =min 0;2^+<y g +x« + ^ + K,2^ 9 + l-2/s;2<^+<y,+x,+ ^+4«,26| fl ;2^ 1 ,+2* 9 +2x g +l+2/c,2^ fl ) , (F4) 



The DSE for the quark-gluon vertex in a skeleton expansion was given in fig. [8l When hard scales are present that do 
not scale to zero there arc different regions of the loop integration that can be IR sensitive. In the uniform case there 
are two possibilities corresponding to a soft or hard loop momentum respectively. In the soft-particle limits different 
propagators in the graphs can involve both soft and hard momenta and there is one additional IR sensitive region 
of the loop integral given by a narrow window of hard modes. This corresponds to a second, not equivalent way to 
route the large momentum through the loop. When assessing the counting of the vertex, each of these regions could 
dominate and determine the scaling of the corresponding Greens function on the left hand side. Employing the IR 
exponents from the Yang-Mills sector in table |T] yields the equations for the IR-exponents of the quark-gluon vertex 
in the different limits 

1 ..„.__„.. 1 
2 

5f g = min (o; 25* + \ + K, 28\ g + S q + Xq + 3 - 2k, 1 - 2k; 26* + \ + Ak, 25 q qg + S q + Xq + 2, 0; 

28* g + l+2 K ,28^ g +28 q +2 Xq +2,0) , (F5) 
8\ g = min (0; ^+^+^+^+2, 8^+8^+2, 8\ g ; 8^ g + 8 q qg + 5 q + Xq + l + 2 K , ^+i| fl +l+2«, 8\ g ; 

28l g +8 q + Xq +2,S u qg +8f g +8 q + Xq +l+2K,8l g ) . (F6) 

where semicolons separate the contributions from the first four graphs in fig. [H] and commas those from different 
regions of the same loop integration. Since the ghost loop contribution in the first diagram in fig. [8] is IR divergent it 
is immediately clear, and reflected in the power counting above, that the second (non-Abclian) graph is sublcading. 

Although the remaining system of equations for the quark sector looks extremely involved, it is largely simplified by 
powerful constraints arising from it. In particular, the non-linear terms on the right hand sides of the vertex equations 
yield the constraints 

% g + 8 q + X q + \+K>0 , S$, + \+K>0 > 8% g + 6 q + Xq + 2>0, (F7) 

Using these constraints in the above equations obtained from the DSEs as well as the known possible range for k [23| 
they precisely cancel the non-linear contributions from the Abclian graph in the uniform and the soft-gluon limit as 
well as all non-trivial contributions in the soft-quark limit and the propagator equations. Dropping the distinction 
between contributions from the different graphs, the remaining system reads 

8 qv - min (25 qv + 1 , 2S qs ) = min (0, 8^ qgv , 8%) , (F8) 
5 qs - min (28 qv + 1, 28 qs ) = min (0, 8%) , (F9) 

1 
2 



5 u qg =min 0,28 q l g +8 q + Xq +-+K,25 qg + l-2 K ) , (F10) 



8f g =min \28f g + - + K,28 qg +5 q + Xq +Z-2K,28l g +28 q +2 Xq +2,l-2Kj , (Fll) 
8* a =wm{Q,S* B ) . (F12) 

To solve the system let us start with the equation for the soft-quark limit of the quark-gluon vertex. From a purely 
mathematical point of view, the second element in the minimum function yields a condition which is trivially fulfilled 
for any <J| < that is compatible with all constraints from the other equations. However, since the perturbative 
vertex does not contain a soft-quark divergence and there is no dynamical contribution that could induce a singularity 
we conclude that 8% g — is the only possible physical solution. 

Next consider the equations for the propagators. Due to the trivial vertex in the soft-quark limit they reduce to 

8 qv = min (28 qv + 1, 28 qs ) , 8 qs = min (2S qv + 1, 28 qs ) (F13) 
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It is easy to see that this system only has the trivial solution 8 q 



0. Correspondingly the self-consistent 



solution has the same structure as the bare propagator with a mass term that suppresses IR propagation. This mass 
is naturally expected to be enhanced by spontaneous symmetry breaking. 

As can be seen form eqs. (|F10|) and (|F11[) the vertex equation in the soft-gluon limit is totally decoupled from the 
equation in the uniform limit. Therefore, it is possible to solve the soft-gluon equation first which is considerably 
simplified by the trivial solution for the propagators and the vertex in the soft-quark limit. Due to the absence of 
propagators in the equation there are no restrictions from the Dirac structure in this case. Thereby the exponents for 
the scalar and vector part are equal and it is not necessary to decompose the equation into different tensor components. 
The remaining equation reads 



^=min(25£ + i+s,l-2« 



and has the two solutions 



= 1 - 2k 



(F14) 



(F15) 



The first one is a strong IR divergence that is self-consistently enhanced by the quark dynamics and which pre- 
cisely saturates the above bound, whereas the second one is the mild singularity fed in from the corresponding soft 
singularities in the gauge sector. 

As has been noted before in the main text, due to the different scaling of the two components of the propagators, 
the remaining uniform vertex equation has to be decomposed into two equations for scalar and vector components. 
This decomposition is performed for the remaining IR leading diagram in figs. 1301 and [3T1 




FIG. 30: The vector part in the decomposition of the IR-leading order of the skeleton expansion of the quark-gluon vertex. 
The open blob represents the ghost triangle in fig. [S]that represents the IR-leading contribution to the 3-gluon vertex. 




FIG. 31: The scalar part in the decomposition of the IR-leading order of the skeleton expansion of the quark-gluon vertex. 



Correspondingly, the remaining system for the components of the vertex in the uniform limit is given by 
5 u qgv = min (^25 u qgv +6 qv + 1 +k, 25% gs + 6 qv + 1+k, 5 u qgv + 5 u qgs + 5 qs + ^ + K,l-2n 
<5" = min ( 25" +S qs +-+K,2S" +5 qs + -+K,5 a l QV +6" +*„ + l+/s, 1-2/s 



(F16) 
(F17) 



Aside from the last terms in these equations arising from the mild soft gluon divergence of the 3-gluon vertex, this 
is the system that has been discussed in detail in section fill Dl and it is not affected by the soft singularities. As 
in the case of the soft-gluon singularity discussed above, these latter terms in principle also allow for the possibility 
of a weaker singularity S qg = — 1/2 — n that is not self-consistently enhanced but only induced from the gauge sector. 
However, our numerical results presented in sec. IIVI clearly favor the strong divergences as a solution of the full DSE 
system. 

So far we have distinguished between the scalar and vector components of quark correlation functions since this 
was crucial to obtain the correct IR exponents that properly take into account the Dirac structure. This distinction is 
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surely also important to assess the chiral properties of the above solution and we recall once more that the scalar part 
is crucial. Aside from this, however, now that the solutions have been found these are just different parts of the same 
Green's function. Thus analogously to the power counting in the gauge sector where only the IR exponents of the 
leading tensor structure(s) are considered, the solutions for the IR-leading behavior of the quark sector of quenched 
QCD is given in Table En 



8, 




8 q 







-3 - k V 1-2/t 





-3 - k V 1-2k 



TABLE VII: The IR-exponents for the leading Green's functions of the possible IR fixed points of the quark sector of quenched 
QCD. The scalar Dirac parts are of leading order in the IR. 



APPENDIX G: DIRAC ALGEBRA 

Skipping for now the color index to concentrate on the Dirac space indices, the vertex in eq. (|C2[) can be decomposed 
with the help of Lorentz invariance and the sign convention in eq. (|A2[) as 



1P1 + r 2P £ + r 37 M + r 4 fa P >i + r 5 M + r e M + r 7 fap% + 

rs faY' + r 9 fol^ + TlQ fa faPi + Tn pi JP2P2 + Ti2 fa t>2ln 



(Gl) 



Our computer codes will calculate the vertex as a Dirac matrix each of whose elements is in turn an integral over 
momentum space. Once this matrix is obtained we need to project out the various structures in eq. (|G1[) . Now we 
briefly describe the algorithm that achieves this. First choose (without loss of generality) a coordinate frame in which 
Pi = (pii 0,0,0) and P2 = (j>2 cos%, 0, 0,p2 sin%). Then pick a vector R of unit norm R-R = 1 along the axis 1,2, 4 
alternatively and numerically compute R-T. The possible structures appearing in this projected vertex are now: 



i?-r NA = rii + r 2 fi + Ts ^i + r 4 i> 2 + r 5 fa /R + r 6 h fi + r 7 fa fa + r 8 fa t> 2 fi , 



(G2) 



where R-T^a takes different values after projection and subtraction in the following three step algorithm, like in eq. 
(|G6|G7j) below, and the various Tj are functions of the (Euclidean) scalars (pl 2 ,p2 2 ,pl-p2). 
They can easily be obtained by tracing the vertex R-T with various projectors: 
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with 



A = p\p\ - (pi -P2? + 2 P i -p 2 R-piR-P2 - p\{R-P2) 2 - pI(R-pi) 2 
In case the vector R is in the pi, P2 plane, the system degenerates. Then the expansion 

R f = Tii + r 3 fa + r 4 fa + r 7 fa fa 

suffices, and the components Ti can be read off 



(G3) 



(G4) 
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(G5) 
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If one chooses first a vector along the 2 (or 3) axis, and employs the full projection above, the combinations I^, 
r 8 , r 9 , r 12 can be separated. Then choosing R along the fourth axis picks up in addition the p 2 contributions, and 
projecting now 

r t = r 4 - r 37 4 - r 8 fa^ - r 9 fa^ - r 12 fa i> 27 4 (G6) 

we can reconstruct T 2 , T$, T-j, T\\. Once these are known we can calculate along the first axis to pick up the p 1 ^ 
contributions extracted from 

r t = r 1 - v 2 p\ - T37 1 - r 5 fa P \ - r 7 fa v \ - r 8 faj 1 (G7) 
-r 9 hi 1 - r n fa fap\ - r 12 fa faj 1 

which finally allows to calculate the remaining T\, I^, Tg, Tip. 



1. Degenerate points pi — P2 



The kinematical limit p\ = p 2 := p, which corresponds to the infrared limit of the soft gluon section, cf. IA 21 
requires special attention. Here A defined in eq. (|G3[) vanishes and the linear system becomes ill-defined: this can be 
tracked to a reduction of the linear space spanned by the vertex, now 
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(G8) 
(G9) 



(G10) 



To obtain the various projected functions the algorithm is in all similar to the general case. With p along the 1 axis, 
choose R along the 2 axis (R-T = T 2 ) and solve for the barred T by means of eq. (|G10j) . then equating 



-psy?n j^ s y m 

1 q 1 9 



1 8 — 1 5 



3 — 1 2 

Next choose R along the p (first in our kinematics) axis and employ 



i?-r = r 1 -i7"y 



1 8 



(Gil) 



after solving 



FT = i^(^-r) rr = ^Tr(M-r) 



(G12) 



the last needed terms can be extracted 
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(G13) 



2. Transformation to the Ball-Chiu basis. 



If wc denote collectively the various Li and Ti by Hi, and the tensors in Dirac space with one index [i that appear 
in eq. (jGlj) above by Hi then, since the vertex accepts a decomposition in both basis, we have 

r = hiHi = Tjli . (G14) 

We wish to obtain the new coefficient functions hi (split into the four WT functions Ai ... A4 and the eight transverse 
functions T\ ... t 8 ) in terms of our calculated r, i C (1, 12). This is given in terms of a linear transformation (not 
unitary since neither basis is orthonormal) 
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Observing cq. (|G14[) it is obvious that 



Hi 



but the matrix A { . is simple to obtain by directly expanding cq. (|A2[) and passing 7^ to the right to compare with 
cq. (IGip . The resulting matrix has two separate blocks, one connecting the terms with an odd number of 7 matrices 
between themselves, another with the even terms. Inverting it wc immediately find 
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expansion in the naive basis eq. 



with A = (j>2 — Pi)- We have to keep in mind that there is a singular kinematical limit. When p\ 
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p the vertex 
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and since all the transverse tensors vanish Tj(pl 
naive basis in eq. (|G8|) reduces to 



)y + (r 4 + r 5 + r 6 + r 7 ) ^ + (r 

= the decomposition in the basis of cq 
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H) referred to the 
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APPENDIX H: DIAGRAMMATIC CONSTRUCTION OF THE QUARK-QUARK SCATTERING 

KERNEL 



While one does not need to be particularly careful when constructing a truncation of the heavy-quark-heavy-quark 
scattering kernel, the chiral Ward identities need to be respected if one is to extend such model to light quarks. 
This implies an intimate relation between the scattering kernel and the gap equation and forbids the naive skeleton 
expansion used in the previous subsection. 

A systematic recipe to construct a scattering kernel that is consistent with a given quark propagator, respecting the 
chiral Ward identities, has been given by Munczek [7^. The basic idea is to take functional derivatives of the quark 
self-energy with respect to the quark propagator. This procedure can be iterated for any type of semi-pcrturbative 
or diagrammatic construction of the quark-gluon vertex by cutting all additional quark lines in the contributing 
diagrams. 

In the most general case, i.e. with an arbitrary vertex, Munczek's construction would entail cutting the selfenergy 
diagrams 
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This would naturally lead to a Dyson-Schwingcr equation for the connected quark-quark scattering kernel 




that couples the four-point quark scattering kernel to a hi gher five-point function. This has been recently and 
independently stressed by Matevosyan, Tandy and Thomas [4l| . However, the confining properties of this kernel 
are not apparent. Indeed if we analyze the scaling of the first term on the right hand side as a function of the 
gluon momentum we obtain kernel{k 2 ) ~ (fc 2 ) _1 / 2_K (fc 2 ) 2K_1 ~ (fc 2 ) K-3 / 2 which is not singular enough to repre- 
sent a confining potential. If, however, we would have two dressed quark-gluon vertices in this diagram we obtain 
kernel(k 2 ) ~ (fc 2 ) _1_2K (A: 2 ) 2K_1 ~ (k 2 )~ 2 , leading to a linear rising potential as discussed in subsection IV Bl 

In order to make this behavior of the kernel explicit we propose the following scheme: (1) we employ a skeleton 
expansion for the (partially) 2PI quark-gluon scattering kernel to obtain a selfconsistcnt equation for the quark-gluon 
vertex which reproduces our results from the main part of this paper, (2) we employ a symmetric formulation of 
the quark self-energy to obtain Munczek's quark-quark scattering kernel with appropriate propagator cuts. This 
symmetric formulation will lead to disagreement with perturbation theory, but this can be systematically controlled. 
Then in step (3) the higher n-point-functions appearing in Munczek's construction can be obtained from a closed 
system of equations by cutting the system in step (1). 

Let us give all these equations in diagrammatic form. Step (1) is materialized by giving the exact equation for the 
quark-gluon vertex in t-channel formulation, fig[2J the BS equation for the connected kernel appearing there, and a 
model for the 2PI scattering kernel in terms of the quark-gluon vertex: 





'mr w\ /w \ 




(H2) 



These three equations form a closed system whose input are the dressed quark propagator, address below, and 
the dressed gluon propagator and three-gluon vertex. In the quenched approximation considered in this work these 
functions do not depend on the quark propagator. 

For step (2) let us express the self energy in a more symmetric form. By substituting the quark-gluon vertex DSE 
in t-channel form, figdl into the self energy one obtains 



(H3) 



We now cast this self-energy in a form that is amenable to systematic improvement. 



n 



(H4) 
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This equation is easily obtained by substituting the DSE for the quark-gluon four-point function, second of eqs. (|H2[) . 
into eq. (|H3p above. 

Resumming the series in (|H4[) using J2^Lo nx ™ = x ( t^x) • wcau d bring us back to the asymmetric formulation of 

the self-energy. In the present form eq. (|H4|) is useful since it explicitly contains the confining part of the interaction 
given by the diagram with two dressed quark-gluon vertices. A truncation of cq. (|H4[) at a given n serves as a practical 
model for a confining kernel which can be worked upon using Munczek's construction. The obvious approximation is 
to keep the first term with two dressed vertices and discard all terms in the sum. The resulting self-energy disagrees 
with perturbation theory at the two-loop level. But this is no worse than the often employed approximation of keeping 
the self-energy formally exact but approximating the quark-gluon vertex at one-loop in the skeleton expansion as in 
the third part of eq. (|H2[) . Moreover, the inclusion of terms with a higher n iterations of the 2PI kernel systematically 
improves the perturbative properties of eq. (|H4[) . Note that also the 2PI kernel was approximated at first order in the 
skeleton expansion, so it also disagrees with perturbation theory at one loop (two- loops in the self-energy) . This also 
needs to be systematically improved with the help of the skeleton expansion to finally match perturbation theory at 
any given order. 

Cutting the self-energy equation then yields an equation for the quark-quark scattering kernel: 




(H5) 

where we have used the derivative of the vertex with respect to the propagator yielding a five-point function. To 
obtain an equation for this five point function we take a derivative of eq. (|H2[) with respect to the quark propagator 
and yielding 




(H6) 

It is thus apparent that we have achieved a truncation that is both chirally symmetric and confining. The quark 
propagator can be obtained from the first term in equation (|H4[) with quark and gluon vertex given by simultaneously 
solving the two-equation system (|H2|) . The matching quark scattering kernel is then obtained from equation (|H5P 
with the attending five-point function given by the two equations in (|H6j) . One has two systems of two equations that 
need to be solved sequentially. 

It is also apparent how the truncation can be matched to perturbation theory to any given order in a systematic 
way by expanding a geometric series of 2PI quark-gluon scattering kernels, and expanding the 2PI kernel itself. The 
number of terms as usual grows very rapidly (but also in perturbation theory). This is a hard task, and it is far from 
our intention to develop this program in the near future, but at least we have found an algorithm for the wanted 
construction that may allow statements of principle. 

Alternatively, and especially for heavy quark applications, one is not so interested in preserving chiral symmetry 
explicitly, since it is already badly broken by the large current quark mass. Then a simple skeleton truncation of the 
kernel may represent a suitable approximation. 
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APPENDIX I: COUPLING TO PHOTONS 

In this appendix let us comment on the quark-photon vertex. Hard electromagnetic probes (Deep Inelastic Scattering 
off the nucleon, jets in electron-positron collisions) and final states (Drcll-Yan processes) are the tell-tale experiments 
about the existence of quarks. Moreover, electromagnetic form factors also provide a glimpse of hadron structure at 
intermediate energies. It is therefore worth paying a moment's attention to the quark-photon coupling to see what 
effect if any does the infrared counting have on it. 

The Bcthc-Salpetcr equation for the quark-photon coupling to flavor i, that for Nf active flavors, and in term of 
the j —> i (possibly flavor changing) scattering kernels reads 




The first term is irrelevant for the purpose of infrared counting. The second displays the quark-quark scattering 
kernels, that have an infrared dimension of (p 2 ) -2 This is compensated by the integration measure d p, so that 
the infrared counting becomes trivial. Thus the equation is consistent with an infrared finite solution. This is also 
what one would expect from the Abelian Ward-i dentity and vector meson dominance. In addition, past and recent 
phenomenological studies of various form factors [llli 1 11 2L seem to do well without the need for an infrared 

enhanced vertex. 

A further question is whether any experiment so far is actually sensitive to the infrared divergences. We cannot 
think of an observable that really constrains this. For example, if one thinks of time-like pair production or jets at 
e~e + colliders, the diagram that would show our enhanced quark-gluon vertex is 




One sees immediately that, by the very definition of a jet, the soft-gluon kinematic point is cut out of the data. 

It does appear that on e ca nnot constrain the infrared counting easily from experiment. A first attempt in this 
direction is given in ref. [115| . but much more work is necessary. Similar considerations apply to the coupling of 
quarks to other external currents, such as axial currents necessary for the weak interactions of hadrons. 

APPENDIX J: THE GROSS- YENNIE ARGUMENT AND ITS FAILURE IN QCD 
1. Crossed ladders simplify ko pole analysis in QED 

In the following we explain why the familiar Gross- Yennie construction of a potential for photon exchange in QED 
cannot be applied to gluon exchange within QCD. 

In an equal-time description, all physical particles are put on-shell and energy is not conserved in the interaction 
vertices. Therefore we start by considering the first skeleton diagram of our (heavy) quark-quark scattering kernel in 
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the center of momentum (CoM) frame, with total energy W — P° as customary, 



(f + q°,q) (f + q ',q') 



;? ~q° -q) (f -q°',-q') 



(Jl) 



and apply our IR counting. In Euclidean space, if all scales were simultaneously sent to zero one obtains 
(PjD 2 ( _1 / 2 ~ K )~ 1+2K = -pr i.e. a strong IR kernel in terms of the Euclidean momentum Pe- With + q ) (and all 



others) fixed at its on-shell value by a Dirac delta function (see also below), all other singularities in the +?o) plane 
cancel and one can perform a Wick rotation. The statement is then that all particles going on-shell simultaneously, 
the skeleton diagram would be proportional to the inverse fourth-power of the off-shell momentum P. However, to 
write down a Hamiltonian description we only put the four quarks on-shell, the gluon rung giving rise to a potential 
interaction. Therefore the condition that the gluon also goes on-shell for our vertex becoming strongly enhanced, 

(9°-g°') 2 -(q-q') 2 -0 (J2) 

becomes, taking into account the mass-shell condition for the four quark legs that has been imposed beforehand, 

(^ - q f q 2 = M 2 (33) 

(and three others), 

(£ 9 -i^) 2 -(q-q') 2 ^o (J4) 



which therefore leads to the statement 



that is, a strong confining potential with string tension named o\ (for the first term in the skeleton expansion) and 
Fourier transform to coordinate space linear in the relative coordinate 

V(r) = <Ti r . (36) 

Here o\ includes the Casimir color factor 4/3. 

The string tension can be extracted from our kernel by contracting all its indices with appropriate spinors. One 
can of course not trust the Born approximation since the kernel is strong, and therefore an iteration mandatory. 

Upon iteration of the kernel (following diagram on the left) , the quarks in the intermediate state are now integrated 
over the off-mass shell region 



q k q' q k q' 




P-q P-k P-q' p _ q (p+k-q-q') p _ q , 

This equation is strictly valid only in QED, with fermions interacting via photon exchange. There, it has been shown 
by Gross and Yennie how the addition of the crossed-rung diagram leads to the correct one-body limit (the second 
particle moves in a static potential created by the first particle whose mass is large). To see this, assume that the 
thicker bottom line in the diagrams above corresponds to a heavy particle of mass M — > oo. Then the two propagators 
of the heavy particle in both diagrams can be combined. In the static limit, the physical intuition is that the total 
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energy becomes the energy of the lighter particle plus the mass of the heavy particle that can absorb recoil momentum 
without absorbing energy, P = e + M , with e = q a = q' . For the heavy particle we will also ignore the negative 
energy part of the propagator, suppressed by a further 1/M. With this notation, the difference of the two diagrams 
is given by a common factor times the difference of the internal fcrmion propagators in the heavy mass limit, i.e. 



1 



1 



M 2 — (P — k) 2 — ie M 2 -{P + k-q- q'f 



(37) 
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2M \M -P + k -ie ' M - P - k + q + q' - ie 
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2M \ko — e — ie ko — e + ie 



2-7ri(5(fco — e) 
2M 



Thus, if one of the particle's masses is sent to infinity, the intermediate k variables are fixed to their on-shell value 
e as intuitively expected. The Wick rotation is then possible, and one can take the large M limit for the second 
particle if so wished. Gr oss and Yennie have given the induction step leading to the cancellation to all orders of 
perturbation theory [l!6| . by summing all ladders and crossed ladders. The argument is valid in the presence of 
vertex and gluon-propagator dressing functions, if the naive vertex mass-counting of perturbation theory holds, and 
thus their proof is essentially complete for QED in the static (automatically quenched) limit. 

In QCD, however, there appear extra color factors in each of the diagrams, such that the simple cancellation 
mechanism of Gross and Yennie cannot work. We are aware of no extension of the argument to include color factors 
as well as QCD gluon radiation. The need for including all crossed ladders makes difficult to show the cancellation 
explicitly when one rewrites the perturbation series as an integral equation, 




as the canceling poles are entangled between the explicit propagators and the 2PI kernel. Therefore, 
the full Bethe-Salpeter equation is known that consistently takes the one-body limit correctly. 



(J8) 

no truncation of 



2. The Gross equation 

Given that the Salpctcr equation (see section IVII Bl in the main part of this paper) needs an inconsistent approx- 
imation of the energy kg thus introducing an uncontrollable systematic error, Gross devised an equation inspired by 
the Gross- Yennie cancellation mechanism. The idea is to replace the bona-fide propagators S q Sq from the field theory 
by an ad-hoc Gross function, containing a Dirac delta that puts the quarks on their mass-shell, and subsequently 
evaluating the kernel at the pole thus forced. The Gross function in the center of mass frame for undressed fcrmions 
is simple and natural, 

G(k p) 27Td({s ^ - ^ rl)/2) - 2 ^ (fco) <m 

1 ' ' Sg 1 + S, 1 2P°(AI 2 + k 2 - P 2 /4) ' {Jy> 

However, in the presence of momentum-dependent dressing functions, one obtains a quite tedious two-particle prop- 
agator, where we use the shorthands A± = A((P°/2 ± fc ) 2 — k 2 ) (sim. for B), 

((A + ^7V+ - B+) x ((- /;+ ^7°M- - B_)6(k° - k°) (J10) 
[(2(P + 2k )(B + B' + + ((P°/2 + k ) 2 - k 2 )A+A' + ) - 2(P - 2k )(B.B'_ + ((P°/2 - fc ) 2 - k 2 )A.A'_) 

+ (Po + 2fcoK - (Po - 2k Q )A 2 _) (((P°/2 + fc ) 2 - k 2 )^ + B 2 + )] - 1 

where the on-shell pole k° is now located by solving the implicit equation 

B\-B 2 _ + ((P°/2 + k°) - k 2 )A 2 + - ((P°/2 - k°) - k 2 )A 2 _ = (Jll) 
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(a numeric fixed-point method or Newton iteration is necessary since the B + , etc. functions depend on k°). The 
resulting local potential that would go into eq. fTB")) is again formally the same as with the Born or Salpeter con- 
structions, except for the different pole position, possibly causing a shift in the gluon momentum and obscuring the 
confining soft-gluon divergence, and in addition the resulting non-local kinetic energy for the equation is markedly 
more difficult. 

It is difficult to see how a confining potential arises here except in the limit A — ► 1, B — ► M where fco — * 
reproducing essentially the same result than the one of section IVII Bl 
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